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SYNOPSIS 
of the 

DISSERTATION ON 

TRAN SPORT PHENOMENA IN A CLOUD OP SOLID PARTICLES, 

DROPS OR BUBBLES 

Submitted in Partial Fulfilment of 
the requirements for the degree of 

MASTER OP TECHNOLOGY 
in 

CHEMICAL ENGINEERING 

SATISH KUMAR 

Department of Chemical Engineering 

Indian Institute of Technology 
KANPUR 

The complications in transport phenomena of dispersed 
or particulate systems due to non-spherical particle sizes and 
a distribution in drop or particles sizes have been pointed 
out. A cell model has been used to obtain the creeping flow 
velocity profile around and drag on an oblate (or prolate) 
spheroidal particle located in a cloud of identical oblate 
(or prolate) spheroidal particles all of which are oriented 
with their minor (or major)axes in the mean flow direction. 

The predicted settling velocities agree satisfactorily with 
experimental data for spheres when the spheroidal particle 
eccentricity tends to unity corresponding to a sphere. The 
mass transport rate to a single active oblate (or prolate) 

• spheroidal particle surrounded by a cloud of inactive identical 



oblate (or prolate) spheroidal particles through -which a fluid 

is having creeping motion has also been obtained for large 

Schmidt numbers. Predicted Sherwood numbers agree quite well 

with the experimental values of Sherwood numbers for an active 

sphere in a cloud of inactive spheres. Using Sirkar’s [17] 

velocity field in creeping flow around a test drop (or bubble) 

placed in a random size-distributed cloud of drops (or bubbles), 

a solution of the convective diffusion equation is also obtained 

for high Schmidt numbers. The expected value of the quantity 
— 1/2 

Sh Pe ' for an active drop in a cloud of inactive uniform 
sized drops was found to be substantially greater than that 
predicted byWaslo andGal-Qr [17] except at high volume concen- 
tration of drops. When the dispersed phase hold up fraction 
tends to 2/3, the present solution breaks down. If one replaces 
the size distributed population of drops (or bubbles) by a 
uniform sized population of drops (or bubbles) of size b.^, the 
volume surface mean radius of dispersion, the magnitude of 
possible error in the mass transfer coefficient was found to be 
less than 6 per cent in all possible cases for the creeping 
flow regime with no surfactants being present. 
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Defined by (5.5) (for Chapters 5 and 6) 
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Defined by (5.7) (for Chapters 5 and 6) 

Total interfacial area between dispersed phase and 
continuous phase 

Minor semiaxes of spheroidal particle (for Chapters 1 
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f Defined by (2.49a) (for Chapters 1 to 4) 

F Defined by (3.33) 
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defined by (3.44) 
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z-direction. 
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h Metrical coefficient defined by (2.43) 
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CHAPTER 1 


INTRODUCTION 

In many industrially important operations involving 
multiphase processing systems, mass is transferred between 
a large number of solid particles or drops or bubbles moving 
in a fluid medium and the surrounding fluid. The estimation 
of mass transfer rate in such systems is quite difficult not 
only due to the unavailability of velocity profile around a 
particle or a drop or a bubble in the simplest case but also 
due to the following features: particle (or drop or bubble) 
size distribution; the complex geometry (non-spherical) of 
particles; presence or absence of surfactant impurities in 
systems with drops or bubbles; residence time distribution 
etc. Since most of the above features are going to affect 
the mass transfer rates, a realistic analysis of particulate 
or dispersed phase systems should incorporate these aspects. 
Unfortunately, most available analyses avoid these compli- 
cations by assuming the following: uniform sized particles or 
drops; particles or drops of spherical shape; constant 
residence time etc. No analysis is available on the velocity 
profile and transport rates for shapes other than spherical 
in particulate systems. The aim of this work is to incorporate 
some of these features in the analysis of viscous transport 
in such systems. Specifically, the analysis of fluid dynamics 
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and mass transport in a cloud of oblate (or prolate) spheroids 
has been attempted first since Neale and Nader [l] have 
already shown that a realistic model for unconsolidated porous 
media or sedimented media or packed beds is provided by a 
cloud of oblate (or prolate) spheroids. Chapters 2 to 4 of 
this thesis are concerned with this analysis. 

Regarding the effect of size distribution, Hoelscher 

and G-al-Or [15] had concluded that a size distributed droplet 

cloud may be replaced without much error by a cloud of uniform 

size dronlets of size b_ 0 which is the volume surface mean 

52 

radius of the size distributed cloud. They came to this 
conclusion by utilising Levich's [21] expression for mass 
transfer coefficient k for a single drop in the flow field for 
a size distributed cloud and averaging it over the whole size 
distribution. Later G-al-Or and Waslo [16] determined the 
velocity profile around a drop for a uniform sized cloud of 
drops or bubbles and postulated that these results are valid 
for a size distributed system also on the basis of the 
conclusions of Hoelscher and Gal-Or [15]. Olney [20] on the 
other hand, has shown that the assumption of a uniform drop 
size may lead to significant errors in estimates of mass 
transfer rates. Therefore, on the basis of Tam's [6] solution 
of the flow field for creeping flow around a test sphere placed 
in a size distributed cloud of spheres, Sirkar [19] developed 
the velocity profile around a test drop in a size distributed 
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cloud of drops. Such a velocity profile should be used for 
checking the effect of size distribution on mass transfer 
rates in a cloud of drops by defining suitable methods of 
averaging as opposed to the mass transfer rates obtained by 
■Gal -Or and Waslo [17] using Gal-Or and Waslo[l6] velocity 
profiles based on a uniform sized droplet cloud. 

The second objective of this work therefore is to 
determine the mass transfer rates in size distributed cloud 
of spherical drops or particles using Sirkar [19] and Tam's [6] 
velocity profiles and Chapters 5 and 6 of this thesis have 
been devoted to it. This will enable us to check the often 
used hypothesis 'Replace a size distributed system of spherical 
particles or drops by a cloud of uniform sized spherical 
particles or drops having the size equal to the Sauter mean 
diameter of the size distributed system'. In all of these 
analyses slow viscous flow is assumed so that creeping flow 
equations of hydrodynamics are valid. Further high Peclet 
number is assumed for calculating the mass transport rates 
so that the concentration boundary layer is affected by 
velocity in the region close to the solid or drop surface. 



CHAPTER 2 


EYLR0LYNAMIC3 OS' A CLOUD Oj OBLATE (OR 
PROLATE) SPHEROIDAL PARTICLES 

Analysis of packed beds and sedimented media encoun- 
tered in chemical engineering practice require modelling with 
particles having shape quite different from that of a sphere. 

Such porous media are necessarily anisotropic. Neale and 
Nader [1] have recently proposed that such systems may be 
adequately represented by a homogeneous swarm of aligned 
particles having the shape of oblate or prolate spheroids. 

The biophysical problems of transport of oxygen as well as 
ionic and organic solutes to red blood cells involve transport 
of the dissolved species to a random suspension of RBC 
particles. This random suspension of red blood cells has been 
modelled by Buckles ot al. [2], Chilcote [3] and Colton et al.[4] 
to be a dilute suspension of randomly oriented oblate spheroids. 
Analysis of transport in porous media or concentrated suspensions 
with oblate (or prolate) spheroids has an additional advantage 
in that wide ranging particle shapes between thin circular disks 
(or slender circular cylinders) and spheres can be accommodated. 

Analysis of mass transport rates in creeping flow 
regime between a fluid and a particle cloud requires the 
knowledge of the velocity field. There are several models 
’cell models* or otherwise, which predict the velocity field 
around a sphere in a cloud of spheres. The earliest of the 
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cell models’ is due to Cunningham [22] who assumed that each 
sphere in the cloud would be effectively limited to motion 
within a concentric fluid mass, the boundary of this outer 
fluid envelope being solid with no slip at its surface. 

Cunningham [22] envisaged that this envelope’s geometry could 
be fixed by assuming that either each sphere is within a concen- 
tric spherical fluid mass such that the new spheres overlap 
but only by so much that they just include the whole of the 
fluid or the arrangement of spherical fluid cells around each 
particle is such that this fluid spheres are packed in the 
closest possible arrangement with the fluid in the spaces between 
them being supposedly occupied by fluid at rest. Obviously due 
to a certain amount of fluid being at rest in this latter model 
(model 2), we will obtain a lower limit of the settling velocity 
whereas with the former model (model l) an upper limit of the 
settling velocity is obtained as Cunningham [22] showed, since 
in terms of relative motion either the outer spherical cell is 
moving and the inner particle is stationary or vice versa (both 
being dynamically equivalent). Cunningham [22] preferred the 
overlapping model (model l) which is also chosen here but the 
exact dimensions are to be chosen here on the basis of Hoppel’s 
[5] ’free surface theory 1 since Happel and Brenner {18] have 
pointed out (page 387) that in Cunningham’s model since the 
cell volumes are not mutually exclusive, the size of a represen- 
tative spherical envelope must be fixed by additional empirical 
consideration’. Hoppel [5] proposed a ’free surface’ cell model 
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in which each sphere in the cloud would he effectively limited 
to motion within a concentric mass of fluid such that at this 
outer fluid envelope, frictionless flow exists. The fluid envelope 
of Happel [5] will be much less than the overlapping fluid 
envelope of Cunningham’s [22] model 1. In a spherical particle 
cloud, any sphere chosen at random would be statistically expected 
to be surrounded by a spherically symmetric cluster of spheres. 

If the radius of such a particle cluster is assumed to be 2b, 
then the vector and scalar profiles will exhibit a maximum or 
minimum at a radius b which will be the radius of the free surface. 
The radius of the outer solid envelope in Cunningham’s model 1 
could therefore be 2b, whereas that in Cunningham’s model 2 will 
be approximately b, the free surface radius. Therefore, as shown 
in Figure 1, we assume that, for a cloud of spheroidal particles 
with their polar axes parallel to the mean flow, any oblate (or 
prolate) spheroid taken at random may be assumed to be surrounded 
by a mass of fluid having the shape of the oblate (or prolate) ; 
spheroid and an outer solid surface. Further, this spheroidal fluid 
in envelope in model 1 can be geometrically related to a smaller 
spheroidal fluid envelope which satisfies the porosity require- 
ments of the whole cloud as in Happel’ s [5] free surface model. 

Note that the spheroidal fluid envelope in Cunningham's model 2 
essentially coincides with Bappel’s free surface. Recently Neale 
and Nader [l] have modelled practical porous media by means of a 
cloud of oblate spheroids in which an oblate spheroid is 
surrounded by its associated pore space bounded by a confocal 
spheroidal shell with the pore space being determined by overall ; 
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porosity of the system. Barlier Jeffrey [10] had calculated 
the viscosity of a suspension of ellipsoidal particles by 
assuming the ellipsoid to be surrounded by a large sphere. 

2.1 Flow Field in a Cloud of Oblate Spheroidal. Pa.rticies : 

This chapter is concerned with the problem of deteminin 
the flow field around an oblate spheroidal particle loco. ted 
inside a homogeneous cloud of identical oblate spheroidal 
particles having a void volume fraction of e when a fluid 
flows past the cloud under conditions of creeping flow. The 
mean flow is assumed to be parallel to the polar axis of the 
oblate spheroidal particles ail of whom are assumed to have 
the same orientation. The solution that will be obtained is 
also useable for the estimation of the settling velocity of 
the ensemble of oblate spheroidal pa.rticles if we can assume 
that all particles have identical orientation and that the 
relative velocity of the fluid with respect to the particle 
is equal to the settling velocity of the particle. This last 
one however need not be true in the most general case. 

The mathematical formulation of this problem is most 
appropriately carried out in oblate spheroidal coordinates (p , 

T), 9 ) as shown in Happel and Brenner [18]. The ~r ,^and ^ 
in Happel and Brenner's [18] nomenclature has been replaced 
here by p, t and a respectively. The oblate spheroid shown 
in Eigure 1 is defined by 
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< t > 2 + ( ? 2 + <!> 2 = 1 

in cartesian coordinates and the origin of the coordinate 
system is the instantaneous position of the centre of the 
oblate spheroid. The cartesian coordinates Cx,y,z) are 
related to the oblate spheroidal coordinates p, p and <p by 

z = V (a 2 -b 2 ) Sinh p Cos p (2.0) 

x = V”(a 2 -b 2 ) Cosh p Cos p Cos <p 

y = V"(a - 2 -b 2 ) Cosh p Cos p Sin <p 

where p^O, O^p^ii , 

If c 2 = a 2 -b 2 and w 2 = x 2 +y 2 (2.l) 

then p and p are defined by 

z + i w = C Sinh (p + ip) (2.2) 

If we put 

t = Sinh p 

and a = Cos p (2.3) 

the above transformation leads to the relations 

z = C t a, w = C V(t 2 +1) Y( 1-a 2 ) (2.4) 

where 0 4 ^ ^co , -1 ^ a ^ 1 

It is well known that the solution to the complete 
Havier-Stokes equations for flows in two and three directions 



9 


is very difficult. For flows with, constant density and 
viscosity, however, the differential equations may he somewhat 
simplified by formulation in terms of a stream function 
as shown by Happel and Brenner [18], For creeping motions 
under steady state conditions, the equations of motion may 
be expressed in terns of the stream function as 



(2.5) 


where in oblate spheroidal coordinates the operator 

2 

E is given as 


B 2 = 


G 2 (x 2 + a 2 ) 


[(t 2 +1) 




+(i 


-a 


2 ) 




~b « 2 


] ( 2 . 6 ) 


Happel and Brenner [18] have solved the equation (2.5) 
for slow fluid motion around an oblate spheroid placed alone 
in the flow field by considering a trial solution of the form 


^ = (1 - a 2 ) «(t) (2.7) 

They obtained the following stream function 

\j£, = (1-a 2 ) [“ 2 °i T + 2 °2 T ~ 2 °2 ^ 2+1 ^ Cot ~ lT 

+ C 3 (t 2 +1)] (2.8) 

where C-^, C 2 ,C^ are constants which are to be determined from 
the boundary conditions for the given system. 
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BOUNDARY CONDITIONS : 

As we have considered the oblate spheroidal particles 
to be at rest with the fluid streaming past, the following 
boundary conditions are valid for flow around any particle 
because the normal velocity at surface must be zero: 


at 

P = ^ 

or 

H 

•r\ 

P 


(Particle surface) 


( 2 . 9 ) 


Since the tangential velocity at the particle surface 
is also zero, the following condition is obtained: 


\j & m0 

3 * 


at p = B 


or 


= 0 


at x = % 


( 2 . 10 ) 


( 2 . 11 ) 


3’ a 

(at particle surface) 


The assumptions of our model proposed in the introduction 
of this chapter would indicate that at the solid surface 
bounding the outermost fluid envelope corresponding to p = 
shown in Figure 1, the condition of no slip should be valid. 
Therefore 


t 


2 U ** = 0 


at p = pi 


( 2 . 12 ) 

or at t = t, 
b 

(at solid surface) 

Further since this boundary condition is equivalent 


to 
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^ = |u C 2 (t= + 1) (1 - a 2 ) 

at x = we notice that the boundary conditions (2. 9) » 

(2.1l) and (2.12) suggest a form of the stream function given 
by equation (2.7) so that the stream function (2.8) obtained 
by Eappel and Brenner [18] are applicable to the present 
problem also. Further the three boundary conditions are 
sufficient to evaluate the three unknown constants C^, C ^ 
and in the stream function expression (2.8) and their values 
are as follows: 



C 1 = 2 Z 

(2.13) 


= (i-t 2 ) z 
a 

(2.14) 


Z X 

°3 “ 2 

(2.15) 

where 

U c 2 

Z = 9 

(2.16) 


(x*-l) (T-X) + (t - -r ) a -X 

q dL d V 

(x b + 

(2.17) 


7 = x, + (r 2 +l) Got -1 t, — 2 t t. Cot -1 T: 

D D b a b a 

(2.18) 

and 

X = x - (t 2 -1) Cot -1 T 
a a a 

(2.19) 


The stream function expression (2.8) may then 

be 

expressed as 



H* = | [t a -('c|- 1 ) (Cot -1 T A -Cot -1 i:) 



T * (t 2 + l) 

a ] 

(t 2 +1) 

(2.20) 
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For the Uniting case of only one oblate spheroidal 
particle in the fluid instead of a cloud of oblate spheroidal 
particles, we have 

cd which moans 

G = x - (t 2 -1) Cot -1 x (2.21) 

a a a 

and the stream function expression (2.20) reduces to (S-ee 
Appendix B equations from B. 1 to B.5) 


1 



t 2 -1) _ 1 

(-2- ) - (-& — ) Cot 1 t — . 
t 2 +1 t 2 +1 

(la ) _ (2H) cct -1 T 
t 2 +1 x 2 +l a - J 


( 2 . 22 ) 


which is identical to' the stream function for flow around a 
single oblate spheroidal particle as derived by Happel and 
Brenner [18]. 

2.2 Brag on a Oblate Spheroidal Particle Placed Inside a Cloud 
of Identical Particles: 


The most general expression (which is derived in Appendix 
A) for the force exerted on an axially symmetric solid particle, 
bubble or drop by the fluid in the direction of the mean fluid 
flow is 


F = it fi 


w 3 


- 2 % p 


-JL r E 2 V- , 

■&P L -o J 

“b 5 


TJ 


w 




hh 


& w 




(1.23) 
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The integrals are over the whole surface of the axially 
symmetric solid particle, drop or bubble. For an axially 
symmetric solid particle the second integral term will vanish 
because at the solid surface, 


"bP 


= 0 (boundary condition( 2. 10 ) ) 
particle surface 


Thus for a axially symmetric solid particle, 

force expression (2.23) reduces to 

r _3 vzty 

= 71 11 J w - TF [ ~ ] 6,1 

which is available in Happel and Brenner [18], 


the drag 

(2.24) 


For the sake of convenience, we replace p and r\ in this 
equation by equations (2.3) so that the drag expression (2.24) 
for an oblate spheroidal particle is given by 


1 =3, 


F = 
z 


n ^ -rB 2 W, 

^ J * ~ [ ~ ] 

-1 


da 


w 


at 

t=t (2.25) 

a 

where the operator E 2 , and w are defined by the equations 
(2.6), (2.20) and (2.4) respectively. Now 


E 2 {f} = 


[(t 2 +l) (l-a 2 ) 


~5 2 # 
S ~ 2 


'] (2.26) 


C 2 (T 2 +a 2 ) Sr" O a' 

which on simplification reduces to the following expression: 



C 1 (l-a 2 ) x 
C 2 (v 2 +a 2 ) 


(2.27) 
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Sy (2.4) and (2.27) we get 


2L r sill 

At l -2 J at 
u w ^ 



(•t 2 +a 2 ) (t 2 +1) -2t 2 ( 2't 2 +a 2 +l) 

2 . 2 2 S 

(t 2 +<x 2 ) (t 2 +l ) 2 
a a 


(2.28) 


Substitution of (2.28) into (2.25) leads to the 


following drag expression 

1 


it p. Ci P (l-a 2 ) 

_ - - - - 


-1 


[ (r 2 +a 2 ) (x 2 +l) 
(r 2 +a 2 ) 2 a a 

St 


-2 t 2 (l+2 t 2 + a 2 )]doc 
a a J 


(2.29) 


or 


2% \x C 


i [(-c 2 -!) - 


(1 + X 2 ) 


2 \ 2 


a 


2 % 


Sin (2 Cot' 




( 2 . 30 ) 

where C^ is defined by (2.13) and reduces to the following 
on using (2.16): 

0 1= 2-^ (2.31) 

Here G is defined by (2.17) and is also equal to the following: 


G- = Sinh p a + ( Sinh 2 p -l) [Cot -1 Sinh P b ~Cot 1 Sinh p ] 

Sinh {3, 

5. cosh 2 p (2.32) 

Cosh 2 p b ^ 

Further U may be considered as the ensemble averaged 
velocity of the particle cloud if the cloud were settling 1 :■ 
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in the fluid. For a fixed particle cloud with fluid moving 
around it. 


F 


z 


4 % u U C (1+B) 

a 


(2.33) 


where 

4 

Cosh p 

B = TaiSTp an (2 Oot anh - 3111112 p a 

a (2.34) 

or 

B = 1.0 (2.35) 

Thus 


F 


z 


8nuU0 

G- 


(2.36) 


For the limiting case of only one oblate spheroidal 
particle in the fluid instead of a cloud of oblate spheroidal 
particles we have 

t, — ° o which means 

b r 


G- 

co 


= G 


= x - ( x 2 —1 ) (Cot 1 % ) 
a a a 


|T b -oo 

and the drag force expression (2.36) reduces to 


(2.37) 


F z = 


8 it p U C 


x -(x 2 -l) Cot F x. 


(2.38) 


a ' a a 

which is identical to the drag force expression for a single 
oblate spheroidal particle as derived by Happel and Brenner [18], 
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Expression (2.36) also reduces to the well known 
Stokes' law of resistance for the limiting case of a single 
sphere in a moving fluid in the creeping flow regime as shown 
in Appendix B equations (B.5 and B.6). 

2.3 Ensemble Averaged Velocity; 

We have assumed earlier that the inner oblate spheroidal 
particle is at rest while the outer oblate spheroidal fluid 
envelope moves with an ensemble averaged velocity U. This 
problem is identical to that of an inner oblate spheroidal 
particle moving with velocity U at the center of a stationary 
outer oblate spheroidal fluid envelope. Thus the drag force 
experienced by the single oblate spheroidal particle will be 
same in both the situations, i.e. 

E = -*- — S U ^ (2.36) 

Z V.T 

where the oblate spheroidal particle is moving in a direction 
opposite to that in which the fluid was moving when the particles 
were at rest. 

Let an oblate spheroidal particle moving with velocity 
U q in an infinite medium experience the same force F z experienced 
by the same oblate spheroidal particle when moving with an 
ensemble averaged velocity U in a cloud of particles. Now 
for a single oblate spheroid IV and U o are related by 
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F = 
z 


8 % ji U 0 C 


(2.39) 


T b =CX) 


Using (2.36), (2.37) and (2.39) we can obtain the 
ensemble averaged velocity U in a particle cloud as 


Sinh P a + ( Sinh 2 p & -l) [Cot ^Sinh p^-Cot ^Sinh p ] 
Sinh p, 

Cosh 2 p Q 

U Cosh 2 B -, 

u - b _ 

° Sinh p a - (Sinh 2 p -l) Cot Sinh p & 

(2.40) 


2.4 ' Velocity Profile Hear the Surface of an Oblate Spheroidal 
Particle Placed in a Cloud of Similar Particles ; 

In scalar transport at high Schmidt or Prandtl numbers, 
the scalar transport is confined to a region very close to the 
solid surface over which the fluid is flowing. The solutions 
of convective scalar transport equations require velocity 
profile close to the solid surface. In this section we propose 
to determine the velocity profile near the surface of an oblate 
spheroidal particle placed in a cloud of similar particles. 

The following relations exist between the radial and 

tangential velocity components and the stream function for 

any axisymmetric body with orthogonal curvilinear coordinates 

p, p, 9 (Figure l) ; 

h 

v — 71 




(2.41) 
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and 


T _5fL ^ 

r\ 


(2,42) 


w * P 

* 

For an oblate spheroid h^, h^ and w are defined as 
follows: 

(Happel and Brenner [18] Page 515) 

1 (2.43) 


h = bu = h = 


and 


w 


h 


P 11 CV (Cosh 2 p-Sin 2 p) 


C Cosh (3 Sin p 


(2.44) 


<P 


The 6 coordinate of any point in the fluid close to 

by 

the surface of the oblate spheroid defined^the coordinate p a 
can be expressed as 

p = P a + AP (2 - 45) 

For locations close enough to the surface of oblate 
spheroid, Ap may be considered to be negligibly small compared 


(3 so that one can approximate 

Q, 


and 


Therefore 


hp (P ) $ hp (p a , p) 




(2.46A) 

(2.46B) 


h = h n ^ h 0 = 


P P a C V (Cosh 2 8 -Sin 2 p) 

8 . 


(2.46C) 
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and 

h $ (2.46D) 

< P' C Cosh 8 Sin p 

cl 

The expression (2.42) may he written in the following 


i orm 


T -ifi. 

r) 




w T 

Substitution of h, 


d t 
d p 


w 


(2.47) 

■ and from (2.46C), 

dp 


(2.44), (2.8) and (2.3) in (2.47) lead to the following 
expression 


Y = — [- + C (1-t Cot" 1 t) 

C 2 f(Cosh 2 P a -Sin 2 r]) 2 2 

+ 2 C^ t ] (2.48) 

or 

Q 

V = SULJD. [ i. +n - C 0 Sin p Cot" 1 Sinh p 

* C 2 V(Cosh 2 P a -Sin 2 p) 2 2 2 

+ 2 C^ Sinh p] (2.49) 

Since p=p a function of p may be expanded in 

Taylor’s series around the value of the function at P a as 

f(p) = f(p a ) + Ap f(p a ) + f " (P a ) 

+ (2.49a) 

As AP is very small, we can neglect terms having (Ap) 2 
or higher powers of A.P in comparison with terms containing 
£p thus 

f(p) = f(p ) + Ap f’tpj 

" ' . . ■ ■ ■ CL Cl 


( 2 . 50 ) 
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This leads to 

Sinh p = Sinh (p o + &p) 

= Sinh + £$ Cosh p & (2.5l) 

and 

Sinh p Cot -1 Sinh p = Sinh (8 + Ap) Got -1 Sinh (p + A P) 

3 3 

= Sinh 6 Cot -1 Sinh p„ + &p (Cosh p o Cot -1 Sinh p Q 

3 3 3 3 

- Tanh p ) (2.52) 

By substituting C^, C 2 , from (2.13) , (2.14) and (2.15) 
respectively into (2.49) and using (2.5l) and (2.52), the 
expression (2.49) reduces to 
A Sin n 

7 = 

11 V(Cosh 2 p -Sin 2 p) 

3 

where 

2 U Tanh 8 

A = a (2.54) 

and G- is defined by expression (2.33). 

Now if we define 

y =&P = P - P a (2.55) 

we obtain 

V = ^ 3in Tl y (2.56) 

f(Cosh 2 P a ~Sinh 2 r)) 

The radial velocity component may be determined 
either by equation (2.41) or by substituting (2.56) into the 
continuity equation. The latter approach is easier. 


A 


(2.53) 
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The continuity equation is 


or 


or 


or 


or 




0 


^ A ) + A C 


>p \h h 

^ n 9 




lip h 


) = 0 


(2.57) 


9 


^ ( _!§_ ) A 

75T 1 \ = " * n 'fc * 


Y 

(rHH 

9 P 


Using (2.56), (2.46C) and (2.46D) 

~C*) 7 " 

">T < h - ^ 5 = - 2 ° 2 A Cosh P a isin 2 n )y 

o K n ? 

u 


y 

( hHH = - 2C 2 A Cosh P a (sin 2 r,)y 
U 9 


Y, 


fi- 


ll h 

ri <j, 


C 2 A Cosh ^Sin 2 rjjy 2 + 


where is a constant of integration and equal to zero because 
at y=0 


Thus 


v e - 0 


Y q = -A*C 2 h * h Cosh p fsin 2 Ti) 
p T) <p H a V 7 


By relations (2.46c) and (2.46D) we can rearrange the 
above to 

A Cost) y 2 


7 p " “ 


V (Cosh 2 p o - Sin 2 n) 
a 


(2.58) 
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By using expression (2.54) in (2.56) and (2.58) the 

radial and tangential velocity components may be expressed 

as follows: 

Badial component: 

2U Tanh {3 ‘Cost} * y 2 
y — ~ S 

^ G-V(Cosh 2 p a - Sin 2 p) (2.59) 

Tangential component : 

2U Tanh 8 * Sinp y 

y cl 

11 &V(Cosh 2 p = - Sin 2 p ) (2.60) 

8 . 

2.5 Expressions for a Cloud of Prolate Spheroidal Particles : 
The prolate spheroid is defined by 

(§) 2 + ( f ) 2 + (|) 2 = 1 ( 2 . 61 ) 

The coordinates p and p are defined by 

(z + i w*) = C Cosh (p+ip) (2.61A) 

where 

C = f(a 2 -b 2 ) 

w*= f(x 2 +y 2 ) 

If we put 

x* = Coshp and a = Cosp (2.6lb) 

the above transformation leads to the relations 

i = C x* a and w* = cV((x* 2 -l)) V(l-a 2 ) (2.61C) 

where 0 if and -1 ^ a ^ 1 
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All expressions which have been derived for an oblate 
spheroidal particle cloud in Sections 2.1, 2.2, 2.3 and 2.4, 
can also be derived for a prolate spheroidal particle cloud 
by methods very similar to those employed for oblate spheroidal 
particles in the above mentioned sections. Instead of consi- 
dering the new problem separately, all the expressions for 
prolate spheroidal cloud can therefore be obtained from the 
analogous results for an oblate spheroidal cloud by a simple 
transformation to imaginary values of parameters and coordinates. 
For example, by replacing if by i x and C by -iC in equation 
(2.20) [18], the stream function for flow around a prolate 
spheroidal particle placed in a cloud of prolate spheroidal 
particles may be expressed as 


* — ' * 2 

■* = - U - — [ T* - (t* 2 + 1) (Cot h -1 t* - Coth -1 v*) 

2 Gr a a a 


where 


T* (l- T* 2 ) 

cL 

(1 - T* 2 )' 


(t* 2 + 1) ( Y* - X*) + (t* - Tf) 2 X* 

P* 8j ■ - ■ - & 2 


( 2 . 62 ) 


(2.63) 


Gosh p. 

p # 1 D -j 

= Sinh 2 Sinh2 P a + (Cosil2 P a +l) ( Goth Cos P a 


Goth -1 Coshp, ) - Cosh p. 


(2.64) 


Y*= x* + (x* 2 -l) Goth - - 1 x* - 2 x* x* Goth -1 t* 
do b a b o 

(2.65) 
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X*= x* - (%* 2 +l) Goth -1 t * (2.66) 

a 3 3 

and 

w*= CV((t* 2 - 1)) *f(l-a 2 ) (2.67) 

For the limiting case of flow around a single prolate 
spheroidal particle in an infinite fluid medium, we have 


CD 

and the expression (2.62) reduces to the result obtained by 
Happel and Brenner [18] (Page 155). 

By a similar transformation, the drag expression (2.56) 
and the expressions for the ensemble averaged velocity (2.40), 
the radial velocity (2.59) and the tangential velocity (2.60) 
for a oblate spheroidal particle cloud may also be transformed 
to the following expressions for a prolate spheroidal particle 
cloud 



* 

U 

IF 

o 



8 x u U* C 
G* 


Cosh B -(l+Cosh 2 B )[Coth“ x 

3 cl 

Cosh 

- Sinh 2 p 

Sinh 2 8 -^ 

Cosh p Q -(l+Cosh 2 p ) Coth" 1 

3 3 


( 2 . 68 ) 

Coshp a -Coth“ 1 Cosp b ] 


a 


Cosh p 


2U*Coth P o Cosp y 2 
G* V" ( Sinh 2 p + Sin 2 t)) 

3 


(2.69) 


and 


( 2 . 70 ) 



26 


V* 

T) 

When 00 ’ ex P ress i° n (2.68) for the drag 

reduces to the result obtained by Happel and Brenner [18] 
(page 155) for drag on a single prolate spheroid. Expression 
(2.68) becomes the Stokers law of resistance for a sphere 
when 



By the proposed model, we have thus been able to 
determine the following quantities describing the viscous 
fluid dynamics in a cloud of oblate (or prolate) spheroids: 

(i)- stream function 

(ii) drag force on an oblate (or prolate) spheroid in a 
cloud 

(iii) Ensemble averaged velocity 

(iv) Velocity profile close to the surface of oblate 
(or prolate) spheroidal particle. 

It must be pointed out however that this is the 
first attempt of this kind for a cloud of oblate spheroidal 
particles. The effects of orientation of the particles 
around the assumed mean position is -unknown. Further the 
difference in the requirements of analysis for a cloud of 
fixed particles or free particles has not been considered 
here [23]. We have considered only cell models in the 


21 f G oth Tsinp ^y 

G*V(Sinh 2 B + Sin 2 p) 

cl 
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derivations presented here as well as in the literature 
reviewed for the objectives of this work. 

2.6 Relationship of Outer Solid Envelope and Free Surface' 
with Void Fraction e ; 

The average void volume fraction e is defined as 
the ratio of the total volume of fluid to the total volume 
of the assemblage. It is implicit in such an assumption 
that the statistical average distance between any two 
neighbouring particles in a cloud of particles will be a 
direct function of the void fraction. As e is increased, 
the average interparticular distance will increase. In 
the model proposed in the introduction of this chapter, 
the normal distance between the particle surface and the 
outer solid envelope of the cell fluid surrounding the 
particle is directly related to the statistical average 
distance between any two neighbouring particles in the 
cloud. If this normal distance is Ygp then one can j 

approximate the normal distance between the free surface and I 
the particle surface by Y^p = Y^/2. Thus 

Y PP = T SP - T sp/' 2 (2 - T2) : 

where Y. . - is the normal distance between ith and jth 1 

13 

spheroidal surfaces with the former enclosing the latter. 

Equations (2.72) with (3*23 A) gives the relation 
between the outer solid envelope, the free surface of Happel[5] 
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and the particle surface in the proposed model. The 
rol&tion may be deduced as 


Sinh p^ - Sinh p f = Sinh p^ - Sinh $ 
or 

Sinh p b = 2 Sinh p f - Sinh p & (2.73) 

To determine p^, we assume following Happel [5] 
that the void fraction for the free surface - particle cell 
is equal to the overall average void fraction in the particle 
assemblage with each free surface effectively isolating 
each particle. 

Thus 


(1-e) = 


Volume of one particle 

Volume within each free surface 


(2.74) 


This leads to 


(i) for an assemblage of oblate spheroidal particles of 
dimension P a (See Neale and Nader [1]). 


Sinh p^ Cosh p^ = 


Sinh p Cosh p 

cl ■ ci 


(2.75A) 


(ii) and for an assemblage of prolate spheroidal particles 
(see Neale and Nader [1]) 

Sinh 2 p Cosh p 

Sinh 2 p f Cosh p f = JJZT) (2.75B) 

The appropriate value of p & is uniquely defined 
in terns of the basic dimensions of the spheroid a (semi 
major axis) and b (semi minor axis), i.e. 
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Cosh p 
r a 


1 

m -c-) 2 ) 

d 


( 2 . 76 ) 


Ihus the equations (2.73), (2.75) and (2.76) provide 
the necessary relationship between p^ and the void fraction 
e of the assemblage of an oblate (or prolate) spheroidal 
particles having a particular value of the eccentricity (b/a) 
of the spheroid. The above dealt with model 1 of Cunningham 
adopted here. However in case of model 2 of Cunningham, the 
free surface of dimension p^ will replace p^ so that for a given 
s, from equation (2.75A) or (2.75B), is known. This p f 
will then replace p^ in all the expressions derived in 
sections 2.1, 2.2, 2.3, 2.4 and 2.5. 


2.7 V elocity Profile Around a Particle at n = 71/2 Location : 

In section 2.4 we developed the expressions for 
linearized radial and tangential velocity components around 
a particle placed in a cloud of identical particles. In fact 
those expressions will predict the corresponding velocities 
without any significant error only near the particle surface. 
In this section we propose to determine the actual velocity 
profile around a particle at the radial location r)=%/2. 

Prom equation (2.4l) 



(2.77) 
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or 



_EL ^ l Q w 
w w 


( 2 . 78 ) 


Differentiating the relation which comes after 
substituting a from (2.3) into (2.4) with respect to p which 
leads to following at x\=tz/2 

I— 1 = 0 (2.79) 

QT) 

using (2.79) into (2,78) 

Vp = 0 at n = te/2 (2.80) 

Similarly from (2.42) for tangential velocity and 
(2.20) for stream function , one can obtain the following 
tangential velocity component at the radial location p=n:/2. 

U 0 - = J [ T a (l~ “?) “ (^ ” ^ ^ cot ~ lT a ~ c ° t_lT )] 

(2.81) 

Figure 2 shows the nature of the velocity profile 
between two particles at r>='n:/2 according to equation (2 80) 
and (2.8l). 



CHAPTER 5 


MASS TRANSFER TO A SINGLE ACTIVE OBLATE SPHEROIDAL PARTICLE 
IH A CLOUD OP INACTIVE IDENTICAL OBLATE SPHEROIDAL PARTICLES 

AT HIGH SCHMIDT NUMBERS 

Ror large Schmidt numbers in the creeping flow 
regime, Pfeffer [11] had analytically determined the 
Sherwood number for a spherical active particle as a function 
of the Peclet number and the fractional void volume in a 
cloud of uniform spherical inactive particles by using the 
velocity profile derived from the ’free surface model' of 
Happel [5]. For a random cloud of size distributed spherical 
particles, Sirkar [12] obtained an analytical relationship 
between the Sherwood number of an active particle, the 
Peclet number, the fractional void volume and the size 
distribution function by approximating the velocity profile 
of Tam [6] close to the solid surface. The problem of mass 
transport to a non-spherical particle in a cloud of non- 
spherical particles has not, however, received enough 
attention, lightfoot [13] on the basis of Stewart’s work [14] 
has indicated the formalism to be adopted for obtaining 
the Sherwood number for mass transfer to a single particle 
of complex geometry placed alone in the flow field at high 
Schmidt numbers. Using the velocity profile around an 
oblate spheroid in a cloud of oblate spheroids, which has 
been derived in the previous chapter, and Lightfoot’ s mass 
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transport analysis for _an isolated particle with saooth 
but complex spheroidal geometry an analytical solution of 
mass transport to an oblate (or prolate) spheroid in a 
cloud of oblate (or prolate) spheroids has been carried 
out in this chapter. 


3.1 Convective Diffusion to an Oblate Spheroidal Particle : 

The equation which is applicable for the convective 
diffusion of species i from or to the surface of the particle 
under conditions of constant density and diffusivity (at const 
constant temperature and pressure) is 


Q'. 

yr+ (-.VC,) 


= D. 


C; + R, 


*->0 

Steady State 


1 V*~i 

L 


No Chemical Reaction 


or 


(v.\7c i ) = D ± 0 i (3.1) 

The equation for convective diffusion of species i 
takes the following form in the coordinate system adopted 
in Chapter 2: 



C, 




V 

JL 

h 0 


C. 

iv*. 


o * 


G. 

i 


~d r 



(3.2) 


A comparison of the magnitude of the various'' terns 
in the equation (3.2) for conditions close to the surface 
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so that {3 = {3„ + <5 whore 5 is the thickness of the diffusion 

Cto 

layer leads to the following (on the basis of similar 
approximations by Levich [21] 


*&A rO 

^ 2b 




Vp. , f\. c. • 

^ 1 Oo 1 

"bh 2 4 b 2 


(3.3) 


and 


1 (\> 




D^ c i q.. A o, 
bP 2 6 2 


(3.4) 


Further 5 <^<^b (which is the semi-minor axis for an oblate 
t spheroid) for the condition of high Schmidt numbers. Hence 


V°i N> V°1 

P 2 C)h 2 

as has been shown by levich [21] for a sphere. 
Thus equation (3.2) becomes 


(3.5) 


1 


[ 




*>> C- 


D, L h ^ p 


7 "S C. 'X 2 C. 

JQ. 1 — i ] _ — L 


(3.6) 


TJ u ■ p ~C> 11 - ^ P : 

The equation (3.6) in dimensionless variables nay 
be written as follows: 

7>S. 


2 t‘ 
o it 


D. h 
x 


U r v « 1 = __ 


(3.7) 



35 


where 


% 


h 


0 . - 0 . 

1 10 

C. - 0. 

X 10 


: h Q = h 
P r\ 

/* _ 

P “ U 


V** = 

r\ 


V 

—I 

u 


(3.8) 
(2.43) 

(3.9) 


By (2.55), (2.59) and (2.60), the following relations 
nay he used; 

~C) 

= 


"dp 


(3.10) 


Y P “ 


u 


**_ _ 


2 Tanh p o Cost) y 2 

c ' =-K Cosri*y 2 


Of (Cosh 2 p ? -Sin 2 n) 


(3.11) 


** v 2 Tanh p Sint] y 

Y = = = K Sinn y (3.12) 

Tj U - - ~ • 1 


&V"( Co sh 2 p - Sin 2 n ) 

ci 


where 


K = 


2 Tanh p 


a 


(3.12A) 


&V" ( Co sh 2 p a ~Sin 2 n ) 

By substituting (3.10), (3.1l), (3.12) and (2.46A) into 
equation (3.7), we obtain the following 


Re- Sch ) 2 ][-E Cosn y 2 A 


K Sin n 7 


^ y = by 2 


(3.13) 


Define 


% = % (a) 


(3.14) 
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and a = y/6(q) 


(3.15) 


where 6 is the dinensionless local boundary layer thickness 
[21]. Rewriting the equation (3.13) in terns of a we 
obtain: 


Re-Sch*f(l-(§^'^) 2 ) [ -K Cost] 6 3 - K Sinn 6 2 ^ ] 


,2 

do 


p 2 _ 

l, TC 

d o : 


( 3 . 16 ) 


Here we choose 6(q) such that 


K Cosri 6 3 + E Sinq 6 ; 


Ee-Sohlf(l-(§ftp) 2 ) 

Or 

(3.17) 


holds good, so that 


— + 3a 2 |5 = 0 

da 2 dcr 


(3.18) 


We have thus reduced the partial differential equation 
(3.13) to an ODE (3.18). The boundary conditions for this 
equation are: 


(A) At y=0, n;=0 (At spheroidal particle surface) 

or at o = 0, it=0 (3.19) 

(B) At y = p -{3 co, it=l (at outer spheroidal 

D cl 

solid envelope) 


As o -5^- co , % = 1 

Equation (3,18) can be integrated to obtain 


or 


(3.20) 
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u 

= j* Exp (-t 3 ) clt + K 2 


Where EL^ and are constants of integration and t 
is a dunny variable. Application of the boundary conditions 
yields the following concentration profile: 


% _ _A 10 _ ___A 

^ico io ] 5 (4/ 3) 


J° exp(-t 3 )dt (3.21) 


The diffusional flux to any point on the surface of 


the oblate spheroidal particle is 

■&c. 

J ■ h ( ^t } T=o 

or pC. ■J) v 

J - D i > y =o 


(3.22) 


(3.23) 


Where Y is the distance fron the surface along a 
line of constant p and is a dinensional quantity, while y 
is nondinensionalised Y. The nomal distance Y fron any 
point (p , p , <p) on an oblate spheroidal particle along 
a line of constant p and <p to the corresponding point (p»p,cp) 
on a confocal oblate spheroidal surface is 


Y = C (Sinhp - Sinhp ) 

3 . 


(3.23)A 


which can be reduced to the following by using (2.51) ana 
(2.55) for p being close to p 


Y = C Cosh p y 

cl 


(3.23B) 


^5y^Y= 0 -C Coshp, 


(3.24) 
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The differentiation of C- with respect to y from 


equation (3.2l) leads to 
' (C i~ G io ^/boN 

n(4/3TS« 1 y7=o 


- JL 

"b 7 


«W°io> 1 


1M} XU X / rz ^ tz \ 

y=0 ; • — - n (w) sw 3 * 5 ' 

The diffusional flux expression (3.23) may be reduced 


to the following by using (3.24) and (3.25): 


J = 


D . ( C . — C . ) 

1 ICO iO' 


c (4/3) Co shj3 6(p) 

ci 


(3.26) 


Substitution of the value of K from (3.12A) in (3.17) 
leads to the following partial differential equation for 
<$(t)) : 


n " , ' v x f\ *z Cosh. 6 

; si “ n 'H + 


Re Sch 


or 


I)" 1 


(6 3 Sinh 3 p) = §■ x 


Cosh 2 p G- Sin 2 p 


Sinh p 


a 


After integrating it, one gets 

, , ■ q Cosh 2 p n 

6 3 sin 3 p = ^ x — -x- - — -~ 


or 


^ 3 m-n vs 3 


1 


Sinh P a Re* Sch 
Cosh 2 p, 


Re Sch 


R 

Sin 2 p dp 


sun p = 4 - ginh p - Re Sch l .. 2 

ct 


x a . r « _ sinia-j 

x R<= Snh L 11 ? J 


(3.27) 


Hie 6 variation with p here is similar to that in 
Levich[2l] for a sphere. At p=m, 6 tends to infinity producing 
zero mass flux as in Levich [21]. This will only reduce the 3 
value marginally. 

Substitution of 6 from (3.27) into 'the expression 
(3.26) reduces the diffusional flux expression to 
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P i^ C ico G ±^ alm \ 4 81 ^a Re-Sch- ,1/3 

1 rtL q z z G-l 

T~ ! (4/3) *C Coshp (n- ~~ 2jl ) Gosh P a 

3. <L 

(3.28) 


The Sherwood number (based on the semi-major axis of 
the oblate spheroid) is a dimensionless estimate of the total 
diffusional flux to the surface of the oblate spheroidal 
particle. It can be estimated by integrating (3.28) over the 
particle surface as 


where 

or 


Sh = 

dS = 


C Cosh p 


a 


D.(C. -C- 

1 100 1 


dn "dm 
h. -h 



J dS 
dS 


r) <? 


dS = C 2 Cosh p • Sint] V"(Cosh 2 p -Sin 2 r)) dp *dcp 

St cl * 


(3.29) 


Thus 


Sh 


2% n 

C Coshp,, | J Sinn \ r (Cosh 2 p -Sin 2 p) dp*dm 

a -J w a 

o o 

2 % % 

Di( G ico -G i) X J SlnTi ^^ Cosh2 P a _Sin2 'n) dr ) d <P 

° c (3.30) 


By substitution of the value of J from (3.28) into (3.30) 
and on simplification we get. 


Sh = (|) ly ^ 


' " I 


1 r Re • Sch- ,1/3 

1473) L G J 


p(p a ) 


(3.31) 
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0.85423 [Sa^sStji/S F(p ) 


( 3 . 32 ) 


where 


sf [ ^ ,1/; V(0osh 2 p a -Sin^) d„ 


F(f3 a ) 


Sinhp, 


o M Sinhe Coshp +1 

Sinh 2 '' 5 p a 0o S h 2 /^ a [0othp a+ ~2~- WcSifag ^ 


( 3 . 33 ) 


and G is defined by (2.32). The factor F(j3 a ) has been 
evaluated numerically using a computer program given in 
Appendix D where the integral was evaluated by Trapezium 


rule. 


Since it is preferable to define the Sherwood number 


and Reynolds number with respect to the diameter for a 
sphere, similarly the Sherwood and Reynolds numbers are now 
defined with respect to the major axis of the oblate spheroidal 
particle. Thus 


0.85423 x 2 r Re ♦ Sch -,1/3 Wr n 

TTf 7 % L g J 


= 1.356 [44] 1/3 F(B ) 


( 3 . 34 ) 


Generally the Reynolds number in packed as well 
as distended beds is defined with respect to the superficial 


velocity U where 
s 

U s = O.e 


( 3 . 35 ) 
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where e is the void volume fraction. 

U *a 

further define Pe = -|j — (3.36) 

1 

By (3.34), (3.35) and (3.36) we get 


Sh = 1.356 [ P({S a ) (3.37) 

Por a single oblate spheroidal particle in the 
flow field, the expression (3.37) reduces to the following 

Sh = 1.356 (-|^-) 1/3 F(s ) (3.38) 

co 

where, as in relation (2.37), 6 is defined by 

CO 

G- co = Sinh p a - (Sinh 2 p a -l) Cot -1 Sinhp a (3.39) 
Limiting Case: Single Sphere : 

. Id ■ , 

When p — co (corresponding to Vi) and 

3- Sr 

c 0, P(p ) reduces to the following 

£L 


p(p a ) = 


2 simT 75 ^ J ()) _ aa£a)i/3 




J 


Sin 2 n 


dp 


P(P a ) 


0.80515 

Sinh ^' 5 p 


(3-. 40) 


a 


(as shown in Appendix B) 


— ,4 

3 Sinh p 

a* 

(as shown in Appendix B) 


and 


G 


(3.41) 
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By (3.40), and (3.41) the expression (3.37) for mass 
transfer to a single sphere in a fluid becomes 

Sh = 1.356 x 0.80515 (| Pe) 1 / 5 (3.42A) 

or 

Sh = 0.992 (Pe) 1 / 5 (3.42B) 

which is Levich's classical solution for a single sphere[2l]. 

The limiting solution for a clout of spheres is 
difficult to obtain analytically. One may numerically 
determine the Sherwood number for an eccentricity value of 
— = 0.99 for obtaining results for a cloud of spheres. 

3.2 Convective Diffusion to a Prolate Spheroidal Particle ; 

All the required expressions for a cloud of prolate 
spheroidal particles can be obtained from those for the 
oblate spheroidal cloud by simple transformations pointed 
out in Section 2.5. 

For high Peclet and low Reynolds numbers, the Sherwood 
number (based on the major axis of the prolate spheroid) 
for a prolate spheroidal particle placed in a cloud of 
inactive identical prolate spheroida.1 particles is therefore 
given by 

Sh = 1.356 ' [-^ l/3 P* (p o ) 

£ Gr 


(3.43) 
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where 



and 


Sin 2 r) 


71 

J f(Sinh 2 p + Sin 2 r)) dr] 

o J ( „ - Sifi-n )i/3 


Oosh. 2 ^e Sinh 2 ^e [Tanhe + Cosh E 2 
r a a a a 

Cot -1 Sinh p a ) (3.44) 


Cosh p, •, 

G = 2 Sinh p + (1 + Cosh 2 e ) [Cotla Cosh p 

Sinh 

-Coth ^ Cosh p^] - Cosh p a 


(2.64-). 



CHAPTER 4 


RESULTS AND CONCLUSIONS FOR a CLOUD OF SPHEROIDAL 

PARTICLES 

In this chapter, >e will compare the results of the 
present work on transport in a cloud of spheroidal particles 
with other analytical or experimental results available 
in literature. In the preceding chapter, an analytical 
expression for estimating the diffusional flux to an active 
spheroidal particle in a cloud of inactive spheroidal particles 
was developed. The same expression may be used also for 
predicting the diffusional flux to an active spheroidal 
particle in a cloud of active spheroidal particles. Unfortu- 
nately no work is available either in terms of hydrodynamics 
around an oblate (or prolate) spheroidal particle in a cloud 
of identical particles or in terms of the mass transfer rate 
to or from the oblate (or prolate) spheroidal particle in a 
cloud of identical particles (active or inactive). As the 
present work covers a x-d.de range of particle shapes lying 
between the limits of thin circular disks and slender circular 
cylinders without being valid for these two limiting cases 
(see comments given later), one can compare existing results 
for a spherical particle cloud with the results of the 
present work when eccentricity of the spheroidal particles 
tends towards unity, since a value of eccentricity of 0.99 
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means that the spheroidal particles are almost of spherical 
shape. Great difficulty was encountered in obtaining an 
analytical solution for the particle shape of sphere when 
b=a because of some indeterminate terms and therefore the 
comparison is being carried out with spheroids of eccentricity 
b/a=0. 99. 

In Figure 3, the nondimensional ensemble average 
velocity U/U q as per model 1 has been plotted against the 
void fraction e according to equation (2.40) for various 
spheroidal shapes having different eccentricities. Hiappel [5] 
on the basis of the ’free surface 1 model has derived an 
analytical expression predicting U /U Q for a spherical particle 
cloud. Richardson and Zaki [25] have offered an expression 
on the basis of their experimental data for U /U ^ for semi- 
mentation of Monosized spheres. The results of these two works 
have also been plotted on the same Figure 3, One notices 
that Happel’s [5] results are considerably lower than the 
experimental results of Richardson and Zaki [25]. On the other 
hand, for an eccentricity equal to 0.99, the results of the 
present work as per model 1 are quite close to the experimental 
results of Richardson and Zaki [25] for both low as well as 
high particle concentrations, whereas model 2 results are much 
closer to those of Happel[5] as expected. 

In Figure 4 the quantity normally used in mass transfer 
in packed beds and dispersions Sh.Pe"*^^ has been plotted 
against the void fraction e for various spheroidal shapes 
having different eecentri cities. A considerable amount of work 




mwmnrm 

ipaisiiii 

■mm 


wmm 


Iter's work [ll] 

abelas et a! [ 28 ] 

tenes and Kramers[ 
son and GfcanKoclis 


work tor, mass | : rc 
analytical -works (1 



48 


exists for the diffusional flux to or from a test particle 

placed in a spherical cloud of other spherical particles. 

Pfeffer [ 11 ], on the basis of Happel’s [5] flow field around a 

sphere in an assemblage of spherical particles, obtained an 

analytical solution for the estimate of the total diffusional 

flux to or from a sphere in an assemblage of spherical 

particles. Pfeffer 1 s [11] results are also plotted on the 

same Pigure 4. For an eccentricity eq ual to 0.99, the results 

from equation ( 3 . 37 ) have also been plotted against the void 

fraction e defined by equation (2.74). It may be observed 

-1/3 

that at all particle concentrations the curve of Si Pe ' 
versus e obtained from the present work for (— )= ,99 quite 

£L 

close to those of Pfeffer 1 s [11] . solution. The solution of 

Sirkar [12] which is based on a different approach to the flow 

pattern in a random swarm of fixed spheres also resembles this 

curve of Sh.Pe - ' 3 '^ versus e closely for high values of e. 

Recently Tardos et al. [ 26 ] (who have derived an analytical 

-l/3 

expression for Sh.Pe ' based on the flow field in a bed of 
uniform spheres proposed by Neale and Nader [27]) found a 
similar tendency. The recent measurements of the overall rate 
of mass transfer to a single sphere in a close packed cubic 
array of inert spheres (e = 0.26) by Karabelas et al. [28] may 
also be compared with the results from the present derivation. 
Wilson and Geahkoplis [29] on the basis of their experiments 
for a dilute bed of active and inactive spheres offered the 
following correlation: 
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Sh = Pe 1 / 5 (4.1) 

e 

It is very clear from Figure 4 that a very good agree- 
ment exists between these experimental results and the results 
of the present work. 

Thus the expressions derived in the present work 
predict scalar and vector transport satisfactorily when the 
spheroidal particle eccentricity (b/a) tends to unity corres- 
ponding to a cloud of spheroidal particles. One may therefore 
expect that the present solution will also predict reasonably 
well the scalar and vector transport for 'a particle cloud 
with particles having an eccentricity other than unity since 
the present solution predicts exactly the hydrodynamics of 
a single oblate (or prolate) spheroidal particle in the flow 
field. For example, it has been shown in Chapter 2, that the 
drag on a flat thin disk as well as that on an elongated rod. 
can be obtained from the present solutions for a single 
particle in the flow field. It nay be possible to compute 
the Kozeny constant K for creeping flow through an assemblage 
of elongated rods with their axes parallel to the mean flow 
and compare these with those of other cell model solutions 
available in P-395 of Happel and Brenner [18]. 

For thecase of mass transfer, it is known, however, 
that there are no creeping flow solutions for Sherwood number 
for flow around a flat thin disk placed perpendicular to 
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flow [l3]* Fo attempt has therefore been made to check the 
validity of the existing solution for this limit- Further 
the effect of the sharp discontinuity in flow directions at 
the disk ends as well as the very tortuous paths of fluid 
flow within an assemblage of disks at high concentrations are 
not known at this time in so far as convective diffusion 
analysis is concerned. Thus the limits of applicability of 
the present solution are not available as the particle concen- 
tration increases and the oblate spheroid tends to a disk. 



CHAPTER 


MASS TRANSFER IN A CLOUD OF DROPS AT HIGH SCHMIDT NUMBERS 


A transfer of mass takes place between the dispersed 
phase and the continuous phase in various chemical engineering 
operations such as distillation, liquid liquid extraction, 
spray drying etc where there exists direct contact between 
the fluids. The study of mass transfer in such systems 
requires the consideration of the effects of drop size distributioi 
presence of surfactant impurities, residence time distribution 
etc. Padmanabhan and Gal-Or [24] and Yaron and Gal-Or [30] 
had determined the total interfacial mass transfer rate in an 
assemblage of drops dispersed in another liquid by probabili- 
stically summing up the mass transport contributions of drops 
in every size range. These two studies were based on the fluid 
mechanics of an assemblage of round drops of uniform diameter i 

I 

as carried out by Gal-Or and Waslo [17]. Recently Sirkar [19] 
developed the velocity profile around a round test drop placed 
in a size distributed cloud of round drops on the basis of 
Tam’s [6] solution for the flow field around a test sphere 
placed in a size distributed cloud of spheres under conditions j 
of creeping flow. Therefore, the objective in this Chapter | 

is to derive an analytical result for the mass transfer 


coefficient to a single active spherical drop inside a random 
homogeneous isotropic cloud of inert spherical drops using 


Sirkar' s [19] velocity profile. This will en$l|l$ US' t<5*|_fp|eck 

' CENU - ,/liRY 
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the often used hypothesis ’Replace a size distributed system 
of spherical particles or drops by a cloud of uniform sized 
spherical particles or drops having the size equal to the 
Sauter mean diameter of the size distributed system 1 . Here 
we restrict ourselves to a cloud of round drops having a 
. uniform residence time. Further, high Peclet numbers and 
low Reynolds numbers are assumed so that existing analytical 
techniques are sufficient to obtain the required solution of 
the convective diffusion equation, lastly, a pure system 
uncontaminated by surfactants is adopted for the sake of 
simplicity. 

5 . 1 Governing Equations for Mass Transport and Its Solution : 


For viscous flow around a spherical particle or drop 


(or bubble) of radius b, if the Schmidt number ( i>/Dj 


is 


large and fluid is in creeping flow, a thin concentration 

boundary layer is formed around the sphere or drop (or bubble) 

from the front stagnation point. Since the Peclet number 
U 2b 

(— — ) under such conditions is very large, the concentration 


distribution for species i in the concentration boundary layer 
may be described by the equation [21]. 

"SC . v. 


a 7*c. 'Vc. 

v *+ D. O 3 

r -gr r -£,0 


( 5 . 1 ) 


1 

At high Peclet numbers, the concentration change is 
completed in a thin region very close to the interface of 
the continuous fluid and the drop. This permits us to use the 
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velocity profile valid in this region. We assume that the 
thickness of the concentration boundary layer, 6, is very 
small compared to the radius of the active drop (or bubble) , 
b. This is consistent with the high Peclet number assumption. 
Following Levich [21] (page 405) the variables (r,&) of the 
spherical polar coordinate system are changed to ( ©) 

where is the stream function of the flow. Further an 
axi symmetric flow is assumed so that there is no $ dependence. 


The flow situation and coordinates are shown in Figure 5. The 
convective diffusion equation (5.l) then reduces to 


c)C. ~rS ^ C. 

‘a* 4 " sin2& 557, [b 


(5.2) 


For the continuous phase flow outside the round test 
drop placed in a size distributed cloud of round drops, 
Sirkar [19] had obtained the following velocity profile: 


v = Cos© [U + — + — (2 - 2e“ ax -2ar e"*“')] (5.3) 


\ o 3 

r a r 


V a = [-U + — - 
© o _3 _2„3 

r a r 


A / -ar,„„ -ar,„ 2 „ 2 „-a r 

(e +ar e +a r e -1/J 


(5.4) 


The quantity a is defined by (See Tam [6] or Sirkar [19] ) 


lia 2 = n(b) 3> r (b) db (5.5) 

where n(b) is the number density distribution of drop sizes 
such that n(b) db is the number density of drops with radii 
between b and b*-db and D r (b) is the coefficient depending or. 
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the drop size b, the statistical properties of the cloud 
and the two fluid viscosities p and p where p refers to the 
drop phase viscosity and p to the continuous phase viscosity. 
To illustrate, the drag force in the jth direction on a drop 
of radius b on the basis of the point force approximation 
leads to 

Dragi = D r (b) U.. (5.5A) 

j dir 

where U . is the unperturbed fluid velocity as seen by the 
drop in the jth direction at the centre of the particle of 
radius b in sense of Tam [6], Sirkar [19] has determined the 
unknowns U q ,A and B in equation (5.4) in terms of the quantity 
*a* : 

U o = + f d+ *£■)] (5.6) 



W = [1+ab+ (l+ab + +-s|4)} (5.10) 

and f = p, (5.11) 

The quantity ab has to be obtained from a solution of 


the following equation 
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co 


2v2 


6 % 


f bn(b) [ 


l+ab+ 


a"b 


2ab 


5 + f ^3 + "jT * 3 


2 V, 2 


a b 


a 3 b 3 ■ 


1 + f(l+ ^j) 


-] x 


(5.11a) 


Since D (b) is givenby 


2 1,2 


6% ^ [ l+ab+ jjj-t f (£ + + a-il- + 


2ab , a 2 b 2 , a 3 b 3 


D r (b) 


•)] 


1 + f (1+ 2f) 


(5.11b) 

The stream function for the continuous phase around 
the round test drop was also determined by Sirkar [19] to be 

U o r 2 Sln6 B Sin 2 © A Sin 2 ©: 


y-- 


+ - — - - [ e -ar ( 1+ar ) -1 ] 

( 5 . 12 ) 

Substituting r=(b+y) in the expression (5.12) for the 
radial coordinate of any point in the concentration boundary 
layer and expanding it in terns of powers of y leads to 


y ■ 


U (l+ab) Sin 2 © 

— rr x [ 2bf y + (3+2f) y 2 ] 

4 [ 1+f ( 1+ ^j)] 


+ Terns with higher powers + 
of y 


(5.13) 


The terns containing y 3 and higher powers of y nay be 
neglected as the Peclet number is assumed to be very high and 
we are concerned with very small values of y close to the drop 
surface. Thus 


, , , U (l+ab) Sin 2 © 

S' =-7EI^T^7 x[2My + (5+2f)y2] 


(5.14) 


db 


and 
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b v, 


0 ■ " r Sin© 




c> f 

U (l+ab) Sin 2 © 


Sin© 


o 


— r — x [2bf + 2(3+2f)y] 
4[l+f(l+ %] 

(5.15) 


Substitution of y in terns of 


U (l+ab) (3+2f) 

*v 9 = - V+ 9 


i+f (i+ ®|) 


fron (5.14) in (5.15) leads to 
) r (z-\jj) (5.16) 


where 


b 2 f U (l+ab) Sin 2 © 
z = ° SbT 


4(3+2f) [l+f(l+ *j)] 

Substituting b Vg_ fron (5.16) into the simplified 
convective diffusion equation (5.2) leads to 


(5.17) 


- D. t> 2 Sln 2 & r[-a 


U (l+ab) (5+2f) 


1!> 


© 


i + f(i+ ^|) 


-1 


.X. 




(5.18) 


A completely general analytical solution for this 
type of P.D,E. is difficult as has been pointed out by Gal -Or 
and Waslo[l7]. For the limiting case of solid particles or 
for fluid dispersed phase with very little internal circulation 
(i.e. f-^0), equation(5.13) becomes , 

C^> c -i 


= -D. b 2 Sin 2 © V" [5U (l+ab )] 




-5 e- - --1 ~ J 

T ( 5 . 18 A) 

This equation is identical to that used by Sirkar [12] 
for obtaining the dimensionless diffusional flux to a single 
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active spherical particle in a random spherical inactive 
particle cloud at high Schmidt numbers. 

For the other limiting case of rapid circulation 
inside the drops, the quantity f has a very large value so 
that and equation (5.18) reduces to 

r) C, 33. U f b 3 Sin 3 © (l+ab) ~^ 2 C, 

-=r - -x - . — — — ^ (5.19) 

C> & 2[1 + f(l + 3)] 

Define ©’ = % - © because the concentration boundary 
layer starts growing from © = % towards 9 = 0 (See Fig. 5) in 
a direction opposite to that of the positive © coordinate. 
Using equation (5.6) and changing the 9 variable to ©' leads to 


C>C, D. U o f (l+ab) b 3 Sin 3 ©* 'Vc. 

1 x s 1 



— x j_ 3 

_ 

2W 

(5.20) 


The boundary conditions are 


(a) 

C ± = 0, at _^> = 0 

(at the interface) 

(5.21) 

(b) 

°i = °10 at y= 0 

(front stagnation point) 



and ©* = 0 


(5.22) 

(c) 

°i =C io 

cd (far from the drop) 

(57 2 - 5 ) 


(corresponding to 

the concentration in the 

incoming fluid) 


For the sane boundary conditions, Levich [21] has solved 

a partial differential equation very similar to the quation 

U f (l+ab) 

(5.20). If in equation (5.20) the quantity — § — xtj is 

d n > 

replaced by U , the total diffxs-ional flux to the drop interface, 

S'' 

1^, was obtained by Levich [21] to be 



59 


h = 8 °io b2 (5.24) 

In the present case, therefore, the total diffusional 
flux to the interface of the test drop can be easily shown to 
be 


_ D U (1+ab) f , 
8 ^ 


(5.25) 


The Sherwood number nay be considered as the dimension- 
less diffusional flux to the total interface of the drop. 
Therefore the Sherwood number defined with respect to the 
diameter of the drop is obtained from (5.25) as 


Sh = 


2 V b 


4 Tib D. C. 

x 1CD 


or 


Sh = 


V(3n) 


f (1+ab) 


x 


2U *b * 
S 1 "T 


D. 


or 


~ f (1+ab) 2U b , 

Sh = [ ^ x — — — x SchP 


f(3%) 

Using (5.6), we get 
2 


f (l+ab) 


-o 


(5.26) 


Sh = 




2 U b , 

h x ( — — ) x gch]* (5.26A) 

1+f (1+ % V 


For a single drop (or bubble) in an infinite fluid 
(i.e. drop concentration, C — > 0 and ab — y 0), the equation 
(5.26A) reduces to 

Sh = 0.653 Re 2 Sch* 

■which, is identical to Levich' s [21] solution. 


(5.26B) 
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5. 2 Interpretation of Mass Transfer Coefficients in a Cloud 
of Drops (or Bubbles) : 

In the preceding section 5.1, the mass transfer 
coefficient for mass transfer to or from an active drop (or 
bubble) moving in the continuous phase has been obtained 
when all other drops (or bubble) are inactive. Wien all 
other drops are active, the effect of concentration wakes 
from previous drops will affect the boundary conditions f or 
the present problem. It is not known how one could solve the 
problem exactly for such a situation in a cloud of drops. 

The existing approaches are similar to that of Gal -Or and 
Hoelscher [15]. These authors had evaluated the total rate of 
mass transfer in the whole vessel, by carrying out an 
integration over the entire surface of the swarm of drops in 
the following manner: 


A _ 

N m = J N a dA 
o 


(5.27) 


T ~ d X, A 

where is the average diffusional flux to the surface of 

a single drop of size b such that dA is the total surface area 

of drops in the size range of b to trt-db. These drops constitute 

a subreactor as proposed by Gal-Or and Re snick [31]. For 

example, for spherical drops in the size range of b to b*-db, 

3 V 

dA = — dC where C is the dispersed phase volume fraction 
and Vj is the total volume of the vessel. Gal-Or and Hoelscher 
[15] had related the total interfacial area A between the 
dispersed and the continuous phases, with a suitable size 
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distribution function n(b) for drops (or bubbles), void 
fraction and the average drop (or bubble) radius. By such 
an analysis, they had evaluated the average nass transfer 
coefficient for the continuous phase with size distribution 
and compared this average mass transfer coefficient with the 
average mass transfer coefficient evaluated on the basis of 
a uniform size of the droplets in the cloud when each drop 
(or bubble) is replaced by a drop (or bubble) having a size 
equal to the volume surface mean radius b^ of the dispersion 
defined by 


co 

J b 3 n(b) db 

_ o 

b -, 0 = 

32 co 

Q S b 2 n(b) db 


(5.27b) 


A similar -approach is being adopted herewith regard 


to the mass transf er c oeff icient derived in the preceding 
section on the basis of Sirkar's [19] velocity profile. However, 
the difference between the two analyses are: Gel-Or and 
Hoelscher [15] used Levich’s [21] equation for mass transfer 
coefficient obtained for a single drop in the flow field where- 
as the mass transport analysis in our case is based on the 
velocity profile through a cloud of size distributed drops. 

The assumption that all the drops (or bubbles) are exposed tc 
the same interfacial as well as bulk concentration loads to 
the following expression for the total mass transfer rate from 
the continuous phase to the interface per unit volume of the 
contactor vessel: 
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I b n(b) db 


(5.28) 


where n(b) is the drop (or bubble) size distribution function 
such that n(b)db is the number density of drops (or bubbles) 
with radius lying between b and (b+db). Substitution of 
from (5.25) into (5.28) leads to 


N T = » C i0 V(f V> J 


® U (l+ab) b 3 . 


n(b) db 


(5.29) 


The total interfacial area per unit volume between the 


dispersed phase and the continuous phase, A^, is, then 


■^■rp — 


' f 


4 Ti b 2 n(b) db 


(5.30) 


The overall average mass transfer coefficient for the 


size distributed bubbles (or drops) may be defined as 


® = c i0 A t 

Using (5.29) and (5.30) leads to 

J? U (l+ab) b 3 t 
1 [-B — u ] n(b) at 


k sd = ' f { m > 


b 2 n(b) db 


(5.31) 


(5.32) 


Substituting U from (5.9) into (5.32) leads to 
^ co _ . c 1 


4Df if' -fig L (1 b 3 n(b) db 


^SD ~ ^ 27 7i jx 


b 2 n(b) db 
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or 


where 


K SD “ R SD 


4 Vf(P'-f) g 

f( — ) 

27 ft p. 



CO 

f(^ 1+ |' --^ ) b 5 ) n(b) *db 

O 

CO 


(5.33) 


(5.34) 


J b 2 n(b) db 
o 

Here the factor ab for a particular drop size b is a function 
of the dispersed phase concentration c the viscosity ratio f 
and the size distribution function (n(b)/N) as given in 
equations (5.11 a,b). 

However, if we replace the size distributed system by 
a system of uniform drops (or bubbles) of size b^ 2 , derived 
from the same size distributed system, the value of 'a* will 
change to ’a 1 . Therefore the average mass transfer coefficient 
in the Gal-Or and Hoelscher [15] sense with such a hypothetical 
cloud will be 




°i0 4 ” ” 


(5.35) 


* 3 : 


Substituting for 1^ with b=b^ 2 from (5.25) into (5.35) and 


using (5.9) leads to 


K, 


TJD 


" ^UD ^ 27 


4 Df( *- ) g 


ft p. 


) 


(5.36) 
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( l+ab-r 2 ) 2 

- VC — ] (5.37) 

V (a b 32 ) 

The quantity 'a' in equation (5.34) is defined by equation 
(5.11a) and (5.11b) for a particular. size distribution function 
n(b)/N. We can assume, for example, that the drop size 
distribution n(b)/l may be approximated by that proposed by 
Bayens [32] for coagulation of hydrosols and subsequently 
modified by Gal -Or and Hoelscher [15]: 

f(b) 4 [ -~ ]^ b 2 exp [-o^b 2 ] (5.38) 

where . 



“o = C 


4 j 2/3 


V* (%) 3 ‘ 

co 


and 


b 3 f(b) db ] 


1/3 


(5.39) 


(5.40) 


o 

The quantity IT which is the total number of drops (or b 
bubbles) per unit volume may be related to the volume concen- 
tration, C, of drops by 


H = 


3 3. 


(5.41) 


where b^ is the volume mean radius defined by (5,40). 

Substitution of N from (5.41) into (5.38) leads to 


n(b) 


30 

^3 


( ~ ) 5/ ^ 2 b 2 exp [ -a Q b 2 ] 


(5.42) 
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For a size distributed cloud of drops (or bubble), 
tiie required value of 'a* nay be obtained from (5.11a) and 
(5.11b) for the nunber density distribution function defined by 

(5.42) . The integration limits of distribution function are 
preferably the maximum and minimum drop (or bubble) size 

in any distribution. For the distribution function (5.42), the 

S, 

cumulative percent volume of drops (or bubbles) with b^5b^ 

covers the range 0.062 per cent to 99.965 per cent. Thus one 
can use these integration limits without any significant error. 

The factor a b^ 2 in equation (5.37) can be obtained 
from equation (5.11a) and reduces to the following as given 
by Sirkar [19] for a uniform size droplet cloud: 

[f(i- f)] + - a B§ 2 [i + f (i - f)] 

-a b 32 [(6+4f) ^ ] - 22. (6 + 4f) = 0 (5.43) 

If one knows the viscosity ratio of two fluids and the 
volume concentration of drops, the factor a b^ 2 for a cloud 
of uniform radius b ^2 can be obtained from this cubic equation 

(5.43) . 

After obtaining the value of *a’ and *a r , it is 
possible to find the- coefficient Rg^ in equation (5.33) and 
in equation (5.36) for different sets of values of f, the 
viscosity ratio of two fluids, C, the volume concentration of 
drops (or bubbles) and given values of b^ 2 . ® 1US » regard 
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to mass transfer in creeping flow one nay compare and 
Rgp to check the often used hypothesis ’Replace a size 
distributed system of spherical particles or drops by a cloud 
of unifom sized spherical particles or drops having the size 
equal to the Sauter mean diameter of the size distributed 
system’ . 

Another way to look at the effect of the size distri- 
bution is to consider the mass transfer coefficient to an 
active drop of size b inside a cloud of inactive drops. 

The cloud of inactive drops nay be considered to be 
either size distributed or having a unifom size equal to b^^ 
of the size distributed system. 

The average mass transfer coefficient, K^, for an 
active drop of radius b, placed in a inactive size distributed 
cloud is obtained from equation (5.25) as 


K SD - 


°10 h 2 


or 


where 


K - R yri tlk (f'-fl-Rl 

•a-n-n «-n V L07 -rr it J 


"SD SD * l 27 


“a. = f 


ir p 


(5.44) 


(5.45) 


Similarly, one can also define the average mass 
transfer coefficient, E^, for an active drop of radius b, 
placed in an inactive cloud of drops of unifom size equal 
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to bj 2 of the size distributed system: 

h 

"hjD - 


C ±0 4u b‘ 


or 

where 


^OD “ R UD ^ f-27 X 

3, 


4 fb(p f -P) J 


1 + a L 0 

vrr •— « »■«*»■ <Djrt» j k 

UD ” 1 L lr /-r 1 a 
¥(ab 52 f ) 

b 32 


7t |1 

32 3 


(5--) 


(5.46) 

(5.47) 


The factors l ab ! and 'a b^p * in equations (5.45) and (5.47) 
respectively may be determined exactly in the same way it was 
determined earlier in this section. Thus one can compare the 
value of Rgj. and for different sets of f, C and different 
values of b and can determine whether a size distributed cloud 
may be replaced by an uniform sized cloud for mass transfer in 
creeping flow. 

A recent experimental work by Tan et. al.[20] on the mass 
transfer coefficient in a packed bed of active spheres involved 
three different diameters of benzoic acid spheres. Such a dis- 
continuous distribution described in Table 1 of Tan et al. [32] 
may be represented by 

n(b) = -^r 5 [x 0> ' 51 6(b-0.5l)+x 0>71 6(b-0.7l)+x 1>15 6(b-1.15)] (5.48) 
4it b~ 

3 

where Xq ^ is fraction of spheres having the radius 0.51 cm 
in the mixture which had spheres of two other radii, 0.71 cm 
and 1.15 cm, and 6(b-0.5l) is Dirac delta function with a value 
of unity when b=G. 51 and 0 when b^0.51. The relative fraction 
of spheres of each size are available in fan et al. [32], Such a 
distribution function will be useful for comparing the predicted 
effect of size distribution with that observed experimentally by 
Tan et al. [32] in the other limit of a cloud of solids particles, 
or almost solid -like drons. 



CHAPTER 6 


RESULTS AND CONCLUSIONS FOR A CLOUD OP 
SIZE DISTRIBUTED DROPS (OR BUBBLES) 

The expression (5.26) derived in the preceding 
chapter should be compared with the existing work on the value 
of the Sherwood number for a. single active drop (or bubble) 
inside a random cloud of inactive drops (or bubbles) for high 
Schmidt numbers and low Reynolds numbers. Unfortunately very 
little work is available for a size distributed cloud of 
drops (or bubbles), in the literature under conditions of no 
surfactant contamination. 

Gal -Or and Waslo [17] had made the analysis for the 

present problem but on the basis of a uniform size of droplets 

— 1/2 

in the cloud. The value of the quantity Sh Pe ' calculated 

from expression (5.26A) for different void volume fractions 

on the basis of uniform size has been compared with the sane 
— 1/2 

quantity Sh Pe ' calculated from Gal-Or and Waslo 1 s [17] 
work, for two different ratios of viscosity of continuous 
phase and viscosity of dispersed phase. It is obvious from 
Table 1 that the estimate of the quantity Sh Pe - "^^ from 
the present work is somewhat less than the value of the 
similar quantity from the expression given in Gal-Or and 
Waslo [17] based on the ’free surface cell* model, especially 
at high volume concentrations of drops (or bubbles) where 
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TABLE I: 

COMPARISON OP PRESENT 

AND GAL-OR’S MODEL E0± 

JCDELJO! 
1 A CLOUD 

ffl WASLO 

OE DROPS 

X 

OF. UNIFORH. SIZE 



Volume concen- 
tration of 

Sh re" 

----- 


dispersed 
phase, C 

f=10 

f=100 


Present Wasio-Gal- 

Present 

Waslo -Gal -Or 


work Or [17] 

work 

[17] 

0.00 

0.6225 0.6225 

0.6497 

0.6497 

0.05 

0.7069 0.7932 

0.7337 

0.8181 

0.10 

0.7373 0.8654 

0.7635 

0.8888 

0.15 

0.7634 0.9288 

0.7890 

0.9507 

0.20 

0.7900 0.9898 

0.8156 

1.0100 

0.30 

0.8415 1.1149 

0.8554 

1.1300 

0.40 

0.8850 1.2540 

0.8946 

1.2700 

0.50 

0.9284 1.4217 

0.9405 

1.4363 

0.60 

0.9969 1.6366 

1.0040 

1.6488 
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in the Unit of the dispersed phase hold up fraction approach- 
ing 2/3, the present solution breaks down. The non-applicability 
of the present work for high value of dispersed phase holdup 
fraction follows fron the point force approximation and the 
assumption of no drop-drop (or bubble-bubble) correlation in 
Sirkar’s [19] derivation cf velocity profile through the 
cloud of size distributed drops (or bubbles). It is appropriate 
at this stage to notice that Sirkar [12] found a sinilar 
tendency when he compared his results for creeping flow mass 
transfer to a single active sphere in a random spherical 
particle cloud with those fron the work of Pfeffer [ll]. 

The expressions (5.33) and (5.36) of the preceding 
chapter predict the overall mass transfer coefficient for the 
size distributed active cloud of bubbles (or drops) and for 
the active cloud of bubbles (or drops)of uniform size ^ 32^32 
respectively. The expressions (5.44) and (5.46) predict the 
mass transfer coefficient for the active test bubble (or drop) 
placed in aa. inactive cloud of size distributed bubbles (or 
drops) and the mass transfer coefficient for the active test 
bubble (or drop) placed in an inactive cloud of bubbles (or 
drops) of uniform size b_ n . Here b „ 0 the Sauter mean radius 
of the size distributed suspension. Thus the predictions 
for a size distributed cloud should be compared with the 
predictions for a cloud of uniform size, b^* so one can 

find the magnitude of possible error if the size-distributed 



71 


cloud is replaced by a uniform sized cloud having the same 
Sauter mean radius. 

Using the size distribution function (5.38) of Gal -Or and 
Hoelscher [15], calculations have been carried out on computer 
(IBM 1401, see Appendix D) for a cloud of drops (or bubbles) for 
three different values of viscosity ratio, f. All the calculations 
are based on the mean volume radius b^=0 . 01 cm. Here it should 
be pointed out that one can also carry out the calcinations 
based on some different value of the mean volume radius b^, but 
the conclusions will not differ as was found out. These results 
are shown in Tables, II, and III where the magnitudes of Rgy, 

Sjjjj, are re P° r ’^ e ^ f° r various cases. In Table II, the 

factors Rg D and Syp In expressions (5.33) and (5.36), respectively, 
have been compared for different viscosity ratios and for 
different dispersed phase volume concentrations. It is obvious 
from this tabd.e that in no ease the possible error is more than 
6 per cent when one replaces the size distributed cloud of 
active bubbles (or drops) by a uniform sized cloud of active 
bubbles (or drops) by a uniform sized cloud of active bubbles 
(or drops) having the same Sauter mean radius b^g. Table III 
Is concerned with an active bubble (or drop) of size, b, placed 
In a cloud of Inactive bubbles (or drops) having the volume 
surface mean radius b^- The factors Rgy and Ryy in expressions 
(5.44) and (5.46), respectively, are reported in Table III. 

Here it should be pointed out that the results will be different 
as the normalized size (b /E^) of the active drop will vary. 



TABLE II: COMPARISON OF Rgp AND R^ FOR AH ACTIVE CLOUD OF DROPS (OR BUBBLES) 
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Calculations shown in Table III are for different viscosity 

ratios and a value of b^ = 0.01 cm. Neither R gD nor are 

expected to be a function of b^. This is supported by calcu- 
lations for other b^ values not shown here. It is clear from 
this table that as the dispersed phase valume concentration 
increases the expected per cent error is also increased if one 
replaces the size distributed cloud surrounding the tost drop 
(or bubble) by a cloud of uniform sized drops (or bubbles) of 


size be 


This behaviour is reasonable because at the lower 


dispersed phase concentrations the effect of disturbances due to 
other drops (or bubbles) is reduced as the inter-drop distances 
has become very large. Therefore the effect of the size distri- 
bution in the other drops (or bubbles on the test drop is muted 
at low volume concentrations. If the test drop has a size other 
than the Sauter mean radius, b^* expected- per cent error 

in R UD changes significantly even if the cloud concentration 
and Sauter mean radius of the cloud remain same. The effect 


seems to be insignificant for test drops smaller than the size 
However, as the test drop size increases relative to b^p* 
the expected percent error in R UP increases significantly. When : 

i 

the viscositjr ratio, f is less than one, the expected per cent ; 
error in Rg^ was found to be maximum at higher dispersed phase 
volume fraction for the test drops (or bubbles) of highest j 

size. On the other hand, for viscosity ratio higher than one the' 
expected error in was found to be maximum with the highest 
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size test drops (or bubbles'' at lower dispersed volume fractions. 

These tables indicate further that there is no case in 
which the expected per cent error in R^p is more than 6 per 
cent when compared to Rgp. Thus the analysis confirms that 
if a size distributed cloud is replaced by a cloud of uniform 
sized drops of size equal to b^ of the size distributed cloud, 
the results would not be much different from that derived by 
using a cloud of size distributed drops (or bubbles). The same 
conclusion was obtained by Gal-Or and Hoelscher [15] on the 
basis of an inadequate expression for mass transfer coefficient 
to a drop as pointed oiit earlier. 

Wo experimental data on drops or bubbles is available to 
check the above theoretical conclusions. However, Tan et al. [32] 
made measurements of mass transfer coefficient in a 5 cm long 
packed bed section of active size distributed spheres. This 
section was preceded by a 10 cm section of inert spheres having 
the same size distribution. Tan et al. [32] concluded that 
their data from Re=2 to 1000 was adequately correlated by "the 
following expression due to Karabelas et al. [28] obtained with 
an inactive cloud of uniform spheres 

4 + [0.29 Re 2y/5 ] 4 d - 

J ( 6 . 1 ) 

if the Sauter mean diameter 2 b^ of the size distributed 
spheres is used in the correlation (6.1) developed with uniform 
spheres. Since 'die depth of the active packed zone is only 5 cm 
and the Sauter mean diameter of each of the size distributions 


e(Sh)(Sc~ 1 / 3 ) 


[1.1 Re 1//3 ] 
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used is around 1.50 cm or more (Table I of Tan et al. [32]), 
we nay assume that the active section is very shallow and that 
the mass transfer analysis for an inactive cloud of spherical 
particles of Sirkar [12] will be applicable here. The relevani 


equation is 


( Sh) ( Sch) 


-1/3 


0.992(l+ab 32 )- 

~~ :i73 — 


(Re)' 


( 6 . 2 ) 


where the particle Reynolds number Re is defined with respect 
to the superficial velocity in a packed bed so that this 

mass transfer asymptote for creeping flow is valid up to a 
particle Re = (U^ 2 b^g/ 1 ))^ 10 [12]. In Table IV, with the 
help of the size distribution iunction, (5.48), values of (Sh)“^^ 
have been calculated for all the four mixtures of Tan et al. [52] 
at two values of Re=2 and 10 for a few Values of void volume 
fraction e(=l-c). The values of e ( Sh) (Sc) from an enlarged 
Figure 2 of Tan et al. [32] for Re=2,10 have been obtained and 
recalculated as ( Sh) ( Sc) in Table IV. One notices, first, 
that up to an e=0.5, the predicted values of (Sh)(Sc)” 1 ^ for 
the four different mixtures are very close to one another. The 
deviations between the predictions for mixtures 1 and 2 on the 
one hand and those for mixtures 3 and 4 on the other hand go 


beyond 10 per cent as the void volume fraction decreases below 
0.45. It is known [12] that solutions based on Tam's [6] 
approach become increasingly inexact in this range of void 
volume fractions. Thus the data of Tan et al. [32] provide a 
limited verification for the theoretical conclusions of this 


study. 
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It may be further observed from Table 17 that the values 
of ( Sh ) ( Sc ) f rom the data of Tan et al. [32] seem to be 
somewhat lower than the calculated values shown in Table 17. 

The deviation increases as the void volume fraction decreases 
to a value of 0.4. This behaviour is also due to the increasing 
inexactness in solutions based on Tam’s [ 6 ] velocity profile 
which blows up when e —> 1 / 3 . However, the deviations for e in 
the range of 0.5 to 0.7 are less than 10 per cent. 

Analytical expression (5.26A) has been obtained for the 
expected Sherwood number as a function of Peclet numbers in 
creeping flow around a round test drop (or bubble) placed in 
an ensemble of round drops (or bubbles) having an arbitrary 

size distribution by utilising the velocity profile derived 

■ ■ * 

by Sirkar [19] on the basis of point force approximation and 

other assumptions of Tam [ 6 ], In the limit of a single drop 

(or bubble) placed in an infinite fluid this expression (5.26A) 

reduces to that of levich [23]. The expected value of the 
- 1/2 

quantity Sh Pe ' in a cloud of uniform sized drops was found 
to be substantially lower than that predicted by Gal -Or and 
¥aslo[17] except at high volume concentration of drops. The latter 
conditions are encountered in the limit of dispersed phase holdup 
fraction approaching 2/3 when the present solution breaks down 
since the assumptions of no drop-drop (or bubble-bubble) 
correlation in Sirkar’ s [19] derivation of velocity profile 
through the cloud of drops (or bubbles) are no longer valid. 



Further this demonstrates the magnitude of possible error in 
the mass transfer coefficient if one replaces the size distri- 
buted population of drops (or bubbles) by a uniform sized 
population of drops (or bubbles) of size b-^. It was shown 
that error is usually small (less than 6 per cent) in all 
possible cases for the creeping flow regime with no surfactants 
being present. 
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SUGG ESTIONS FPL FU T URE WOR E 

1. It would be useful to make experimental measurements 
in the manner of Karabelas et al. [28] of the mass transfer 
coefficient to a single active spheroidal particle in a 
packed bed of similar inactive particles with all of thorn 
having the same orientation. The variation of the eccentricity 
ratio will be particularly desirable. Such measurements 

may also be attempted in an ex-panded bed. High Schmidt number 
conditions are to be maintained and natural convection should 
be eliminated for this viscous transport problem. 

2. It would be worthwhile to develop the equations of 
motion for viscous flow around a test spheroidal particle 
in a cloud of size distributed spheroidal particles with 
the help of the point force approximation of Tam [6] and 
investigate the possibilities of getting its solution since 
its boundary conditions are similar to those for a stationary 
spheroid in a fluid moving with a uniform velocity at infinity. 

3. Since nonspherical bubbles are often encountered 
in industrial practice, it would be useful to model the 
motion for a spheroidal drop and obtain estimates for its 
Sherwood number in viscous flow. 
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APPENDIX A 

BRA& ON A BODY 

Ihe second rank (dyadic) stress tensor, -pj-, for 
Newtonian fluids, is given by 

TT = - e 1 + e (v.b 5 + s ii A (a.i) 

where P = the hydrostatic pressure the fluid would be 
supporting if it were at rest at its local 
density and temperature, T. 

I = unit tensor 
K = bulk or volume viscosity 
p = shear viscosity 

-IS 

A = ra -te of deformation tensor, defined in Eq.(A.3) 

For compressive fluids the pressure p which appears in 
equation (A.l) is not the mean noimal pressure, P m , at any 
point. The mean pressure P m is defined as follows: 

- - f TT 8 =I 

= - J <flxx + TT yy + TT ZZ ) 

(for a Cartesian coordinate system) 

P m = P - K (^7.V) (A. 2) 

The term, in equation (A.l) is based on a 

linear relationship between the viscous portion of the 
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pressure tensor and the rate of deformation (shear) tensor. 
This leads to 


A = \ [ V v + (vv) t ]-|i( ! 7 v) (a. 3 ) 

where (\J V; = transpose of 7. 

For an incompressible fluid, the equations (A.l) 
and (A. 3) can be written in more simplified fonn as follows 


TT = - 5 I + 2 H A 

(A.4) 

A= |[VT + (yf)*] 

(A. 5) 


We shall adopt the orthogonal curvilinear coordinate 
system for an axisymmetric body. The notations for the 
coordinate system are same as those in Chapter 2. 

Thus for an axisymmetric body 


h = h 

P B 


h 


and 


(A. 6 ) 


(A. 7 ) 


9 w 

For our coordinate system (see page 489, Happel and Brenner 
[18]). 

. C) v - c> 


V V- i I -h 

p p L p 


( - )1 

r\ *c)T) h 


~ “ *C! V n ^ ,1 3 

+ i r i h [•= -? ■■ - h v fi -c — (r ) ] 

P t) L ?> P Pen y 
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^ T i w ^ ,1 ,, 

'4 V - h T n Ijp ( h P 


U U h C' 


r ^!a 


+ % \ h[ ^ + s T g ^ ( h » 

"V ■ 

+ I I h [h v c~ 7 ~ ( h ) + h v — (-—-) 
9 9 9 L P DP n^,Ti Y 

> 


(A. 8) 


< Vv)* - ip h blJ ff- +hv ? T^ )] 


Hr 

ir 1 

P P 

h[-%^ 
L bp 

- h v 

P 

Hr 

i i_ 

b v 

i a rJ±-IL 

~ l ~&p 

— h v„ 


P p 

p 



> 


Hr 

i i 

h [~J5- 
L Oil 

+ h v. 


p p 

P 


D P ( h } ^ 

— (-)l 

Dp ( n } J 

blT £)] 




+ i i h [h v ?T~r (— r~) + 

999 P OP *y 


li v 


11 'Dp 


<-£->] 


9 

(A. 9) 

Substitution of C\Jj) and ('\J V)" t from (A. 8) and 
(A. 9) respectively into (A. 5) loads to 

'b 


C?D 


/\ = ip ip h C"^TR' + 11 ( *)] 

- T h “~bb 


■W " n p ^ 


+ Vp 2 [ + " h V P ^ ( h } 


+ i i h[ — — ^ + h v — ~r (r - ) ] 

T] T] L s T) P 3p V J 

+ i i h [h V -rpr ({““)+ h v ~4r (~r~) 

CD cp cp 6 B 6 h T) dT) & 

+ Vp 2 L J„ + S p - % ^ ( h ) - h T fi JnV-l 

The stress vector acting across an element of 


rl> 


(A. 10) 


area whose outer normal is i rt shall, be of the fons 

P 


surface 


85 


TT p =H,ip = (-P I + 2\i&) « i 

The unit tensor, I, is given by 


(A. 11) 


I = i„i r +i i + i i 

P P h T) cp (p 


(A. 12) 


Substitution of A and. I from (A. 10) and (A. 12) 
respectively into (A.ll) and further simplification leads 
to 

_ _ 7)V„ "X v n 

Tip = -ipP + 2ph[ i p -^+ i c- ] 

P P L P O P r) q J 


P 


“h y 


+ i n h*j>^r- (^-) - 
rj p L Qf h ' 


IL (An 

>n 1 h } J 


+ 2p v [i - i ~ — ] 

^ T] L T] >p p r) J 


(A. 13) 


The general properties of any vector leads to the 
following two relations: 


V 


T , ~S T C - , v r 

v r = i n h - v" T + i h — — 1 L - 

P P OP h 


(A. 14) 


md 


"">> v r 


Curl V = i h 2 [“7 r (■#•)] 


9 


L *^P 


,n v h 


V (our i f) = h2[ ^r (Z i 5 - It ( -^ )] 


where V has three components v ,v and v (=0). 

H T1 f 

Using equations (A. 14) and (A. 15), the stress 
vector defined by equation (A. 13) reduces to 


TT p = -ip* + 2(1 V v p + % t 1 [y V)] 

+ 2(1 v [i U.16) 

^ r\ L r\ 



Happel and Brenner [18] have derived the following 
relations for a axi symmetric body 

Vorticity vector = Curl V 

Curl V V = B 2 \jJ (A. 17) 

Finally, the relations (2. 41), (2.42) and (A. 17), 
connecting the normal velocity, tangential velocity and the 
vorticity, respectively, to the stream function, leads to 
the following stress vector: 



The force exerted on the body by the fluid in the 
positive z-direction will be 

- Jrr z • 45 = jTh.Ip is 

S= Surface S 

of the body 

■ JH, « - 43 

s s 

2 % 

= lJtL 2 I 4 n d 9 

o T 

F z = ^jTIpz T a n < A - 19) - 

TT pz -n,s. = -<yy p - 2. <yv> <| 

+ ~ | - ( y v Hf (A ' 20) 
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For relating the coordinates, the following relations 
hold good: 


1 1 

w T) 


and 


~§r\ 


1 P 1„ 

P Z 


(A. 21) 


1 ip = h-=r-7r = - 1 1 

w p 0 P r\ z 


(A. 21 ) 

Substitution ofTT p 2 from (A. 20) in to (A. 19) with 
regard to the equation (A.2l) and (A.22) leads to the following 
drag expression 

p 

, W‘ H ) On + 77 1 


=-^ ( w 2 P) <*, + TfJ 


w =r— 6-n 

7) ^ 


w 

p 


6p 


^-v 1 6,1 - ^ f \ 

-**f [^r ‘iBr + -lf <rr 1 sr, 


(A. 23) 


J T)P L ~bP "bp ~ir) -^n 

Where the integrals are to be taken around a Meridian 
section of the body, the end points of the path of integration 
being the upper and lower points at which the boundary surface 
intersects the axis of symmetry. Since w = 0 at these terminal 
points the first a nd third integrals vanish identically, leaving 

_ P-2 *C>P 


= M ^ “ 2*^ E 2 6 1 ) 


J 


-4np.J^ 




w "^> h 


lw 

"dP 

1>w ~*>h-,^ 


] ^TfT 


(A. 24) 


~0P ~T>P "^fl "W'&P 

Happel and Brenner [18] had developed the relationship between 
the pressure in steady creeping flow and the stream function 


as follows: 
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Jtt (B 'f 5 

w 1 

Substitution into the previous expression yields 


(A. 25) 


= TCp 

•w 


_3 

w r-r- 

A b P 


[ 


] <$ti 


w 


- 2%\x 


h 


'~ ^-£ r "^ V ^_h , ~^_W 




[ DP -?>P + Dti 


"bn 


■] 6p (A. 26) 


This is the most general expression for the force 
exerted on an axisymmetric body when the continuous phase is 
an incompressible fluid. 
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APPENDIX B 

SIMPLIFICATION OF -EXPRESSIONS POR A SINGLE .SPHERE 
MOVING IN AN INFINITE FLUID MEDIUM 


The factor ' G' for a cloud of oblate spheroidal 
particles is defined by equation (2.32) as follows: 

G = Sinh p a + (Sinh 2 p a -l) [Cot -1 sinhp b - 
1 Sinhp, 

Cot Sinh 8 ] - Cosh 2 B (2.32) 

a Cosh 2 j3 b a 

Equations (2.73) and (2.75) indicate that as the void 

fraction e increases, the parameter 'p b * also increases. For 
the limiting case of the void fraction e becoming one, this 
parameter becomes infinite i.e. 

As e 1.0, |3 b — -Ss* 6 - co 

or Sinh p b — Cosh p b — co 

and Cot -1 Sinh P b ■ 0 

Thus, equation (2.32) reduces to 

G = Sinhp - — ( Sinh 2 j3 -l) Cot -1 Sinh 8 (B.l) 

cl Q, 3 , 

The function Cot -1 Sinh p & may be written in the form 
of a series as follows: 


Cot -1 sinh p a = 


1 _ 1 + ___1 

Sinhp a 3 Sin h 3 e 5 Sinh 5 6 

d , d 


7 Sinh P a 


(B.2) 
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Substitution of/ this into the previous expi* ess i°& 
yields the following equation: 

1 1 i 


G = 7Z 


Sinh p a 3 Sinh p a 3 Slnh 3„ 5 

a 


a 


or 


+ 


+ 


5 Sinh p 7 Sinh p 


a 


Gr = 


4 1. 


_ JJ 

3 Sinh 8 15 

a 


i 


Sinh 3 p, 


+ 


12 1 
35 Sinh 5 p a 


(B.3) 


Equation (2.76) indicates that the coordinate (3, 

c 

tends to infinity as the ratio b/a approaches unity. 


P a ^ 00 


as f ~>1 


md 


Sinh p a — ^ co as - ^ 1 . Therefore 


since Sinh.6 is a very large number as ~ — ^ 1, we conclude 

cl £t 


that 


or 


Sinh^a sinil P a as a ~ ^ 


Sinh p i 




Sinh 3 p 4 


as 


a 


(B.4) 


Thus the equation(B. 3 ) can be written with the 
help of (B.4) in the following form 


G = 


4 


3 Sinh 8 

a 


as — 
a 


(B.5) 
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Substitution of G from (B.5) into (2.36) gives the 
drag experienced by a sphere moving in an infinite medium 
as 


9 % p U C Sinh 8 

— C 

4 

3 


as 


a 


i 


C Sinh p^ = r, radius of sphere 


when - = l 

cl 


'Thus 


F = 6 % u U r 
z p 


(B.6) 


which is the well known Stoke' s formula. 

The factor l(p ) in equation (3.37) for a oblate 

EL 

spheroidal particle is defined by (3.33) as 


% 


Sin 2 n 




(t) - 


sl y T j j VS r<Oosli 2 p a -Sin 2 T,) dr, 


Sinh 2/3 p 0Dsh 2/5 p [Doth R + (lhhh-0] 


or 


Cosh 7 


n 

1/3 e. S 


JiaflL 


a 


r(p a ) = 


0 (n- Sa2a)i/3 


3 e v Co sh p a -l J 
(3.33) 

ni-( )2 ] 


a 


Sinh 2//3 6 [Gcth p„ + 

cl cu* 


Sinh p Gosh 0+1 

2 ~ log e ( — )] 


Cosh p -1 
r a 


(B.7) 


As we know that 

1 + 


2 l0g e ( 


Cash p 


a 


■) = 


1 


1 - 


Gosh p 


a 


Cosh p a 3 Cosh 3 p 5 Go sfcTp 

a ,. ir 

(B.8) 


a 
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Substituting the above into the previous expression 
yields the following: 


; IT 

^ X 


Ooth 1/3 e f — 


F(p a ) = 


(n - gg- ) 1/3 


Vn-f ) 2 1 dn 

u± 1 Cosh p/ J at1 


Sinh 1 / 3 8 [ Cothp^ + Tanhp + 


Tanh p 


S, 


Tanh p. 


3 Cosh 2 p 5 CoslTp 

ct 

(B.9) 


+. 


a 


Now P a — ^ oo 
and Sinhp 


as 5 — ^ 1 


a 


Coshp 


a 


co as - — 1 
a 


Tanh p 


rsj 


Coth p ^ 1.0 as - — ^1 


a 


a 


and 


Tanh p 


a 


Tanh 8 

l_Q 

Cosh 2 p 


a 


Thus the equation (b.9) can be written for the 
limiting case of a single sphere in the fluid in the follow- 
ing form: 

n 

f _§iah d „ 




p(p a ) 


p(e a ) 


0 SiL2a)i/3 


2 Sinll 1 / 5 S ; 
1-6103 


2 Sinh 1 / 3 p. 


0.80515 

Sinh^ 3 p, 


(B.10) 


(B.ll) 


a r a 

where the value of the integral in the numerator of (B.10) is 


available on page 85 of Levich [21], 
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APPMDTir c 

“ TTI >DCT 

For lower particle concentrations we may assume that 
all particles are uniformly distributed and are spaced 
according to the well known Body Centered Cubic Structure 
for spheres. For such an assemblage of oblate (or prolate) 
spheroidal particles, we assume that the particles are spaced 
according to the BodyCentered Kectangular Parallelepiped 
(BCBP) unit cell as shown in Figure 6 such that only the free 
surface of each particle touches the free surface of another 

particle, while the particles are not in direct contact with 
each other. 

Here we define 


Dispersed phase volume per unit 

1-e = BCRP cell 

Volume of BCRP unit cell 

V 

1 “ e = ~v (c.i) 

Two unequal arms l g and 3^ of BCRP unit cell will 
he in the same ratio as that of the minor to major axis of 
the oblate (or prolate) spheroidal free surface: 

l g G. Sinh 
3^ ~ C Gosh 


or 



■mi 



■■ 

. 


FI EVATION 


Free surface 

Particle 

surface 
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1 s = h IaEl1 P f 

Volume of BOKP unit coll l s g iven by 


( 0 . 2 ) 


V = l 2 2 . 

m s 

using ( C. 2) 

V = X m Tanil Pf (C.3) 

Distance D 0 between the centers of two neighbouring 

free surface envelopes is equal to two times the distance 
between the centre and the point of contact of a free surface 
envelope. Thus 

D e = 2 Cf(Sinh 2 p f + 0.5) (C.4) 

D e = half of the diagonal of the BCEP cell 

= 1V(1 2 + l 2 + l 2 ) 
sms 

This reduces to the following on substitution of 1 

s 

from (C. 2) 


D e = "§• f(2 + Tanh 2 p f ) 

By equating (C.4) and (C.5) we get 
4 Cf(Sinh 2 p f + 0.5) 


1 = 

m 


f(2 + Tanh 2 • p„) 


(C.5) 


(C.6) 


Substitution of l m from (C.6) into (C.3) leads to 


Sinbh 


V 


p f + 0.5 n3/2 


64 C 3 Tanh p f [ ^.0 p~ 1 


(0.7) 
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Each corner free surface envelope does not belong 
uniquely to one unit cell, but is a part of all eight unit 
cells surrounding it. Therefore only one eight part of 
each corner free surface envelope belongs to one unit cell. 
Each cell also possesses a free surface envelope located 
at its center that is not shared with any other unfit cells. 
Thus BCRP unit cell has two free surface envelope --f one 
contributed by the eight corner free surface envelopes and 
one located at its center. In other words we can say that 
there are two oblate spheroidal particles per BCHP unit 
cell. 


Thus volume of solid phase per BCRP unit cell is 

V = 2 x | n C 3 Sinh p o Cosh 2 6 (C.8) 

Substituting V and 7 from (C.8) and (C.7) into (C.l) 
s 

leads to 


Tanh p,, [ 


Slnh 2 jj f + 0.5 , /2 _ £ 


Sinh 


s a Oo^B 


2 + Tanh p_ 


(C.9) 


Eor prolate spheroidal particles assemblage the 
same relationship may be obtained in the following form: 


Coth 


r Sinh 2 p f + 0.5 i5 / 2 
L 2 + Coth " 24 


Cosh p Sinh 2 



(l - e) 


(C.10) 


Thus the equation (2.73), (2.76) and (C.9(lO)) 
provide the relationship between and void fraction 
of the assemblage of an oblate (or prolate) spheroidal 
particles having a particular value of the eccentricity 
(b/a) of the spheroid. 

At higher particle concentration this BCEP model 
will not predict the correct void fraction because the 
probability of violating the basic assumption of BOEP model 
that the particles are not in direct contact increases as 
the particle concentration increases. 
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APPENDIX E 

SIMPLIFICATION OF DRAG EXPRESSION FOR A CLOUD 
OF ELONGATED RODS 


The prolate spheroid resembles a long thin rod if 
the semimajor axis, a, is much greater than the semiminor 

axis, b, such that 

.# 

0 = V(a 2 -b 2 ) (B.l) 

This approximation leads to 


cl cO 

Gosh P a - q ~ 

, „ b 

Sinh P a = q = 


1.0 

b b 
a ” a 


(E.2) 

(B.3) 


The general expression for the drag on a prolate 
spheroidal particle placed in a cloud of identical particles 

is given by (2.68) 



8 tc \x a U* 
G* 


( 2 . 68 ) 


where 


& * = C ° S --h Sinh p a + (1 + Cosh 2 p a ) [CotlT 1 Cosh P a 

- CotlT 1 Cosh p„] - Cosh p a (2 ' 6 ’ ) 

TJsihg (E.2) and (E.3) elation (2.64) reduces to 

» = Cosh jb ^ _ 2 cotiT 1 Cosh p,, + 21h [|J + 2 121 2 ■ 1 ’° 
Sinh 2 p b a (1.4) 
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and 


Cosh, p b = 2 Cosh. p f - 1.0 
Sinh (3 b = 2 Sinh p f -[J] * 


where 


Sinh^p f 



and Cosh p^ = V(l + Sinh 2 p f ) 

^ (1 + Sinh p. - 7) 

X ct 


(E.5) 

(E.6) 


(E.7) 


(E.8) 


Thus tho drug on a elongated rod placed in a cloud of identica 
particles reduces to 


8 u (i a U* 

■p* — • — — 

z Cosh p, -u „ 

[j]-2 Coth~ x Cosh p, + 21n [^]+21n 

Sinh 2 p b a b b 


2 - 1.0 
(E. 9 ) 


For a elongated rod particle placed in an infinite 
medium tho expression (E.9) reduces to 


F. 


4 it u a U 


In [|] + ln2 - 0.5 (E.10) 

which is identical to that obtained by Happel and Brenner [18.] 



